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SUMMARY  FOR  PROPOSAL  62940-CH-REP 

The  current  funded  project  has  been  very  successful  in  the  scientific  development  in  addition  to 
its  broad  impacts  on  the  research  infrastructure,  minority  STEM  student  research  training  and 
community  services.  The  project  was  performed  based  on  the  originally  proposed  tasks  and 
schedule  timeline  with  minor  modifications  according  to  the  actual  process.  Overall,  several 
important  topics  were  investigated  using  combined  experimental  and  theoretical  modeling 
methodology.  Those  topics  are  of  great  interest  and  importance  to  not  only  the  science  and 
technology  community,  but  also  to  DoD  and  Army  Research  Office.  The  research  results  are 
encouraging  and  they  are  reported  as  Part  A-E  in  this  report  in  a  format  of  a  scientific  publication, 
respectively. 
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PART  A:  MECHANISM  AND  KINETICS  OF  PROTON  TRANSFER  IN 
MOLTEN  CARBONATE 

Al.  Introduction  and  Background 

Solid  oxide  fuel  cells  (SOFCs)  are  of  growing  interest  as  clean  energy  source  due  to  the 
advantages  they  provide,  including  high  efficiency,  low  emission  and  fuel  flexibility  Both 
oxygen-ion  (o-SOFCs)  and  proton  conductor  (p-SOFCs)  based  electrolytes  have  been 
extensively  examined.  The  o-SOFCs  normally  operate  at  high  temperatures  of  700-1000  °C, 
causing  severe  issues  like  thermo-mechanical  compatibility  and  chemical  stability,  and  also 
leading  to  high  fabrication  cost.  For  p-SOFCs,  the  operating  temperature  is  much  lower  (400-600 
“C)  using  the  perovskite  of  BaCeOa  (BC)  and  its  doped  derivatives  as  electrolytes  making 
them  promising  candidates  for  commercial  power  generators.  However,  the  cerates  react  badly 
with  CO2  and  become  dysfunctional.  Despite  many  efforts  to  improve  the  performance  of  BC- 
based  electrolytes,  BaZrOa  (BZ)  and  its  derivatives  are  considered  as  alternatives.  The  BZ- 
based  electrolytes  are  chemically  stable,  but  challenged  by  low  conductivity  at  the  desired 
temperature  range.  To  achieve  the  commercial  feasibility,  durable  electrolytes  with  high 
conductivity  at  400-600  “C  is  necessary. 

Yttrium-doped  barium  zirconate  (BZY)  is  one  of  the  promising  proton  conducting  electrolytes 
for  intermediate  temperatures  (IT)  SOFCs  due  to  its  excellent  chemical  stability  and  good  bulk 
proton  conductivity.  However,  its  refractory  nature  makes  itself  very  difficult  to  sinter  and  hard 
to  achieve  high  density  structure.  More  recently,  Huang  et  al.  reported  that  densification  of 
BZY  can  be  promoted  by  adding  molten  carbonate  (MC)  as  a  sintering  aid.  MC  was  filled  into 
the  pores  of  BZY  first  by  infiltration  method.  Interestingly,  not  only  high  density  structure  of 
BZY  was  observed,  but  also  high  ionic  conductivity  of  0.33  and  0.38  S/cm  at  600  °C  in  3%  H2O- 
air  and  3%  H2O-H2,  respectively.  The  contribution  from  proton  transfer  to  the  ionic  conductivity 
was  found  to  be  about  55%  and  the  MC  phase  had  played  a  significant  role  in  the  proton 
conduction.  The  conduction  mechanism  was  mainly  proposed  by  the  equation  below. 

H2C(,)  +  Vo-  +  (C03)^oo3  =  +  (HC03)* 03  Eqn  (1) 

This  enhanced  conductivity  implies  that  the  MC  phase  in  BZY  can  provide  extra  channels  for 
proton  migration.  Similar  hybrid  systems  of  MC  with  BZY,  BCZY  and  even  oxygen-ion 
conductor  as  TT-SOFC  electrolyte  have  been  previously  reported  by  other  groups  and  largely 
enhanced  cell  performance  was  observed.  Such  composite  electrolytes  are  easy  to  fabricate  with 
low  cost,  which  will  open  a  door  for  developing  novel  electrolytes  for  IT-SOFCs.  In  addition,  it 
is  important  to  point  out  that  the  contribution  from  hydroxide  ion  (OH')  to  the  proton 
conductivity  was  not  considered  here,  while  this  was  reported  to  be  noticeable  in  Ref.  25.  In  fact, 
the  effect  of  water  on  the  proton  conductivity  was  already  examined  independently  by  two 
different  groups  Both  reported  that  the  conductivity  was  increased  by  only  6-8%  when 
water  partial  pressure  changes  from  0  atm  to  0.3  atm  in  5%  H2  Ar  or  N2.  This  increase  was 
explained  due  to  the  reaction  of  water  with  the  surface  defect  (oxygen  vacancy)  and  carbonate  to 
produce  HCOa'  and  OH'  as  proton  carriers.  Therefore,  it  is  fairly  reasonable  to  think  the 
contribution  of  HCO3'  to  the  proton  conductivity  is  much  larger  than  that  of  OH'.  The  role  of 
OH'  is  not  examined  in  the  current  study,  but  should  be  considered  in  the  future. 

In  bulk  perovskite  oxides,  proton  transfer  occurs  through  hopping  between  adjacent  oxygen  ions 
at  normal  lattice  sites  via  a  Grotthuss-type  mechanism  Further,  experimental  and 
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computational  data  indicate  that  the  proton  transfer  involves  two  steps:  rotational  diffusion 
around  an  oxygen  ion  and  transfer  diffusion  toward  a  neighboring  oxide  ion  For  proton 
conduction  in  BZY/MC  as  described  in  Eqn  (1),  the  protons  produced  by  surface  defect  reactions 
were  transferred  to  the  neighboring  carbonate-ions  (COs^')  at  the  BZY/MC  interface  to  form 
HCOs'.  They  will  then  be  transported  inside  of  the  MC  phase.  However,  such  assumption  has  not 
been  verified  in  any  form  yet.  Previously,  we  have  reported  a  static  DFT  study  of  the  proton 
transfer  in  the  crystal  structure  of  lithium  carbonate  The  calculated  energy  barrier  was  0.34  eV 
along  the  direction  of  [001].  However,  the  experiments  in  Ref.  16-17  have  shown  that  the  MC 
phase  only  have  large  contributions  to  the  ionic  conductivity  when  it  is  melt  at  higher 
temperature  than  450  °C.  This  implies  that  the  results  in  Ref  33  are  not  enough  for  understanding 
the  proton  conduction  mechanism  in  molten  carbonate.  In  the  current  study,  two  significant 
improvements  have  been  made.  One  is  to  use  a  disordered  cluster  model  to  represent  the  MC 
phase,  which  is  close  to  the  real  situation  and  should  produce  more  reliable  results.  The  second  is 
to  include  the  effect  of  finite  temperature  using  the  first-principles  molecular  dynamics  (FPMD) 
method. 

A2:  Modeling  Methodology 

In  particular,  the  cluster  of  (Li2C03)8  was  used  to  represent  the  MC  phase  and  those  calculations 
were  performed  in  two  steps.  In  the  first  step,  all  gas  phase  geometries  of  [(Li2C03)8H]‘^  cluster 
were  optimized  at  the  B3LYP  ^"^■^^/6-31G(d,p)  level  using  the  Gaussian  09  program  All 
stationary  points  on  the  potential  energy  surface  (PES)  were  then  verified  by  calculated 
vibrational  frequencies  at  the  same  theoretical  level.  Therefore,  the  reactant  and  product 
structures  are  truly  local  minima  while  each  transition  state  at  the  first  order  saddle  point  is  on 
the  PES.  Eor  such  a  large  molecular  system,  locating  a  transition  state  is  time  consuming  and 
challenging.  To  quickly  achieve  this,  we  performed  a  geometry  scan  first,  then  used  the  structure 
with  the  highest  energy  for  TS  search,  and  finally  verified  the  reaction  pathways  using  the 
intrinsic  reaction  coordinate  (IRC)  calculations.  In  the  second  step,  we  applied  the  VASP  code 
(version  5.3.5)  to  perform  first-principles  molecular  dynamics  (EPMD)  calculations.  The 
optimized  structure  of  [(Ei2C03)8H]'^  from  step  one  is  used  as  the  starting  geometry.  All  EPMD 
simulations  were  conducted  with  the  NVT  ensemble.  PAW-PBE  potentials  were  used  for  H 
(ultrasoft),  Ei  (slpO),  C  (s2p2),  and  O  (s2p4)  with  the  energy  cutoff  of  500  eV.  A  Verlet 
algorithm  was  integrated  with  Newton’s  equations  of  motion  at  a  time  step  of  0.1  fs  for  a  total 
simulation  time  of  3  ps,  i.e.,  30000  time  steps  in  total.  The  frequency  of  the  temperature 
oscillations  was  controlled  by  the  Nose  mass  during  the  simulations.  Additional,  a  1x1x1  k-point 
mesh  was  used  at  the  F -point. 

The  proton  transfer  process  is  local  and  the  reaction  is  usually  at  high  temperature  in  a  molten 
carbonate  salt,  both  of  which  will  bring  detriments  to  the  results  from  static  DFT  calculations  in 
the  gas  phase  using  a  limited  size  of  cluster  model.  For  example,  the  calculated  energy  barrier 
includes  contributions  from  not  only  the  proton  migration  itself,  but  also  the  variation  of  the 
cluster  geometry.  Therefore,  there  are  possible  uncertainties  in  the  calculated  energetics.  To 
verify  and  correct  such  possible  errors,  we  have  re-examined  the  energies  from  local  structures 
extracted  from  the  cluster  and  the  revised  values  showed  large  improvements  and  better 
consistency  with  current  FPMD  results  and  previous  experimental  results  In  addition,  long 
range  interactions  in  the  molecular  system  were  estimated  by  CAM-B3EYP  but  no  significant 
difference  was  observed. 

A3:  Results  and  discussion 
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A3.1:  DFT  results  of  proton  migration  in  the  (Li2C03)8  cluster 
Intra-carbonate  migration  of  proton 

For  the  migration  of  proton  in  MC,  both  intra  and  inter-carbonate  pathways  are  possible  and  they 
will  be  discussed  separately.  As  reported  in  Ref.  33,  the  cluster  of  (Li2C03)8  is  fairly  reliable  to 
represent  the  disordered  molten  carbonate  with  affordable  computing  time.  Fig.  1  shows  the 
optimized  structures  of  reactant,  transition  state  (TS),  product  together  with  the  relative  energy 
for  proton  migration  inside  of  a  carbonate  ion  in  the  cluster  of  (Li2C03)8.  At  the  beginning,  the 
proton  (fF^)  is  connected  to  01  with  the  Ol-H"^  distance  of  1.001  A.  It  then  starts  moving  toward 
02,  leading  to  a  TS.  At  TS,  the  only  imaginary  frequency  is  /1851  cm'\  corresponding  to  the 
shifting  of  from  01  to  02.  From  reactant  to  TS,  the  Ol-H"^  bond  is  stretched  to  1.313  A, 
while  02-H'^  shortened  to  1.276  A.  At  the  same  time,  C-01  is  shortened  from  1.354  A  to  1.297 
A,  while  C-02  elongated  from  1.261  A  to  1.328  A.  It  is  also  noticeable  that  the  Ol-C-02  bond 
angles  is  bent  to  106.6°,  which  brings  01  and  02  closer  and  will  facilitate  the  proton  migration. 
The  Mulliken  population  analysis  shows  a  loss  of  0.02  e  for  01  from  reactant  to  TS,  gain  of  0.05 
e  for  02  and  loss  of  0.04  e  for  carbon.  All  results  indicate  that  the  proton  migration  is  a 
concerted  process  and  involves  2p  orbitals  of  01 -C-02.  Fig.  1(c)  displays  the  main  bonding 
orbital  between  proton  with  01  and  02  at  TS.  It  is  very  clear  that  the  2p  lone  pairs  of  01  and  02 
combine  with  H  Is,  forming  a  bent  Ol-H-02  linkage.  The  cleavage  of  Ol-H"^  and  formation  of 
02-IT^  is  reversible  as  verified  by  the  IRC  calculations. 

The  calculated  energy  barrier  for  proton  migration  from  01  to  02  is  46.8  kcal/mol  and  40.4 
kcal/mol  for  the  reverse  process.  The  energy  difference  between  the  initial  and  final  state  is  6.4 
kcal/mol.  With  long  range  corrections  by  CAM-B3LYP,  the  barrier  is  47.6  kcal/mol  and  39.9 
kcal/mol,  respectively,  showing  negligible  contribution  from  long  range  interactions  in  the 
system.  However,  the  value  is  only  20.5  kcal/mol  in  Ref.  33,  which  was  calculated  based  on  a 
single  molecule  of  HC03'.  To  eliminate  possible  uncertainty  in  the  calculated  energy  barrier  as 
discussed  before,  we  have  extracted  a  small  HC03'  cluster  out  of  the  whole  structure  and  re¬ 
calculated  the  energies.  The  energy  barrier  is  decreased  to  17.2  kcal/mol  and  17.5  kcal/mol  for 
the  reverse  step.  The  energy  change  for  the  migration  process  is  reduced  to  0.3  kcal/mol.  Both 
imply  that  the  contribution  from  geometry  variation  in  the  cluster  other  than  proton  migration 
itself  is  significant  and  the  results  from  such  calculations  should  be  carefully  evaluated  and 
verified. 
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Fig.  1:  (a)  Structures  of  reactant,  transition  state  (TS)  and  product,  the  upper  panel  is  the 
whole  cluster,  while  local  structure  only  for  the  lower  panel;  (b)  Relative  energy  for  proton 
migration  inside  of  a  carbonate  ion  calculated  by  the  whole  cluster  and  local  structure, 
respectively;  (c)  Selected  molecular  orbital  in  TS  for  bonding  between  proton  and 
carbonate.  The  proton,  carbon,  oxygen  and  lithium  atoms  are  shown  as  yellow,  gray,  red, 

o 

and  purple  respectively.  (Distance  in  A  and  relative  energy  in  kcal/mol.) 

Inter-carbonate  migration  of  proton 

Fig.  2  shows  the  optimized  structures  of  reactant,  TS,  product  as  well  as  the  relative  energy  for 
proton  migration  through  an  inter-carbonate  pathway  in  the  cluster  of  (Li2C03)8.  Unlike  the 
intra-carbonate  process,  this  pathway  involves  two  neighboring  carbonates.  As  shown  in  Fig. 
2(a),  a  TS  structure  for  the  migration  of  between  02  and  05  was  located  on  the  PES.  At  TS, 
the  only  imaginary  frequency  is  /387  cm'\  corresponding  to  the  migration  of  from  02  to  05. 
At  the  beginning,  the  02-H'^  distance  is  1.041  A,  similar  to  Ol-H"^  in  Fig.  1(a).  However,  it  is 
only  1.492  A  for  05-H'^,  which  is  0.784  A  shorter  than  02-H'^  in  Fig.  1(a),  implying  that  the 
inter-carbonate  pathway  is  possibly  favored  by  structure.  At  TS,  a  linear  02-H-05  linkage  is 

o  o 

formed  between  two  carbonate  ions  with  the  02-H  and  05-H  distance  of  1.152  A  and  1.281  A, 
respectively.  Fig.  2(c)  shows  large  overlap  between  2p  of  02  and  05,  and  Is  of  H.  The  Mulliken 
population  analysis  indicates  a  gain  of  0.10  e  for  02,  loss  of  0.02  e  for  05,  and  gain  of  0.03  e  for 
Cl  from  reactant  to  TS.  A  similar  trend  is  observed  for  the  reverse  path  from  product  to  TS. 
Although  the  absolute  values  of  charge  are  not  very  reliable,  the  small  variations  of  electron 
density  on  the  atoms  where  the  proton  migration  occurs  indicate  that  this  is  a  synergetic  process 
between  proton  and  the  two  carbonate  ions. 

The  calculated  energy  barrier  is  21.6  kcal/mol  from  reactant  to  TS,  while  only  9.6  kcal/mol  from 
product  back  to  TS.  It  becomes  20.8  kcal/mol  and  8.2  kcal/mol  after  the  long  range  correction  by 
CAM-B3LYP,  respectively.  The  large  energy  difference  between  reactant  and  product  (12.0 
kcal/mol)  implies  the  fact  that  the  cluster  of  product  from  IRC  might  be  a  local  minimum  on  PES, 
i.e.  an  intermediate  structure.  A  similar  phenomenon  was  observed  in  Ref.  33,  where  a  very 
small  activation  energy  of  1.8  kcal/mol  was  reported  for  the  proton  transfer  between  two 
carbonate  ions.  Using  an  extracted  cluster  Ei2C03-H'^-Ei2C03  as  shown  in  Eig.  2(a),  the  energy 
barrier  is  reduced  to  8.0  kcal/mol  and  8.3  kcal/mol,  respectively.  This  again  confirms  that  the  gas 
phase  calculation  may  overestimate  the  activation  energy  for  a  molten  salt  system  using  a  simple 
cluster  model.  Erom  the  current  study,  it  is  found  that  the  calculated  energy  barriers  will  have  to 
be  checked  and  verified  by  other  methods,  which  will  be  discussed  in  more  details.  Comparison 
of  the  calculated  energy  barriers  suggests  that  the  proton  migration  proceed  through  inter- 
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carbonate  rather  than  intra-carbonate  in  the  lithium  molten  carbonate,  which  agrees  with  the 
results  from  the  geometry  and  charge  analysis. 
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Fig.  2:  (a)  Structures  of  reactant,  transition  state  (TS),  and  product,  the  upper  panel  is 
cluster  structure  and  the  lower  panel  is  local  structure;  (b)  Relative  energy  for  proton 
migration  for  inter-carbonate  in  the  (Li2C03)8  cluster  and  local  structure,  respectively;  (c) 
Selected  molecular  orbital  in  TS  for  bonding  between  proton  and  carbonate.  The  proton, 
carbon,  oxygen  and  lithium  atoms  are  shown  as  yellow,  gray,  red,  and  purple  respectively. 
(Distance  in  A;  Relative  energy  in  kcal/mol.) 

A3.2:  FPMD  evidences  of  proton  migration  in  the  (Li2C03)8  cluster 
Trajectories  of  the  proton,  carbon,  and  oxygen  atoms 

To  account  for  the  effects  of  temperature  on  the  proton  migration  in  the  lithium  molten  carbonate, 
we  carried  out  first-principles  molecular  dynamics  calculations.  Fig.  3  shows  a  series  of  the 
snapshots  of  proton  position  at  1300  K  in  a  time  interval  of  0.5  ps  and  the  corresponding 
trajectory  of  proton,  carbon  and  oxygen  ion  in  the  (Li2C03)8  clusters  simulated  at  1300  K  for  3 
ps.  The  proton  position  in  each  snapshot  has  experienced  noticeable  change  with  the  time, 
indicating  a  fast  proton  diffusion  process.  The  trajectory  of  proton  is  largely  dispersed  while  the 
trajectory  of  carbon  is  quite  localized,  indicating  that  the  proton  migration  in  the  (Li2C03)8 
cluster  is  very  fast  with  large  displacement,  while  the  carbon  atoms  only  vibrate  around  their 
original  sites  with  very  limited  displacement.  It  should  also  be  noted  that  the  trajectory  of  the 
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oxygen  atom  is  more  dispersed  than  that  of  carbon,  implying  that  the  oxygen  atoms  can  rotate 
and  adjust  their  positions  during  the  proton  migration  between  the  carbonate  ions  in  the  cluster  of 
(Li2C03)8. 


(b) 


Fig.  3:  (a)  Snap-shots  for  proton  transfer  in  the  cluster  of  (LizCOajs  simulated  at  1300  K  for 
3  ps.  The  yellow  ball  represents  the  proton  in  motion,  (h)  Trajectory  of  proton,  carbon,  and 
oxygen  ion  in  the  cluster  of  (Li2C03)8  simulated  at  1300  K  for  3  ps. 

Mean  square  displacement  and  Arrhenius  plot 

The  mean  square  displacement  (MSD)  is  determined  by  the  ensemble  average: 

MSD{t)  =<dr\t)  +  Equ  (2), 

where  N  is  the  number  of  ions  in  the  system,  n  is  the  number  of  time  origins,  t  is  the  time,  and  to/ 
is  the  initial  time  step  originating  at  time  i  In  the  current  work,  at  least  three  FPMD  runs 
were  conducted  for  one  structure  at  each  temperature  in  order  to  obtain  an  averaged  MSD  for 
better  reliability  and  accuracy.  Fig.  4  shows  the  averaged  mean  square  displacement  for  proton 
migration  in  the  (LiiCOsjg  clusters  simulated  at  the  temperature  of  1100  K,  1200  K,  1300  K  and 
1400  K  for  3  ps  (here  considering  the  melt  point  of  lithium  carbonate  is  about  996  K).  Clearly, 
the  MSD  of  proton  increases  linearly  with  time,  indicating  a  fast  proton  migration  in  the 
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(Li2C03)8  clusters.  Moreover,  the  slope  of  MSD  increases  as  the  temperature  increases, 
indicating  that  the  proton  has  higher  ability  of  diffusion  at  increased  temperature. 

The  diffusion  coefficient  is  calculated  by  fitting  the  slope  of  the  MSD  vs.  time: 

Eqn(3). 

6t 

Then,  the  Arrhenius  plot  can  be  obtained  from  the  diffusion  coefficient  via  Arrhenius 
relationship: 

Z)(r )  =  (r ) exp ( — Eqn  (4). 

kgT 

Herein,  and  are  the  pre-exponential  factor  and  the  diffusion  barrier,  respectively,  is  the 
Boltzmann  constant  and  T  is  the  absolute  temperature.  By  fitting  to  the  Arrhenius  relationship 
(Eqn.  4)  over  the  considered  temperature  range,  the  pre-exponential  factor  and  the  effective 
activation  energy  can  be  obtained.  Eig.  5  shows  the  Arrhenius  plot  for  proton  migration  in  the 
cluster  of  (Ei2C03)8.  Erom  the  slope  of  the  straight  line,  the  activation  energy  is  calculated  to  be 
7.5  kcahmol,  which  is  in  good  agreement  with  the  energy  barrier  of  8.0  kcal/mol  for  the  proton 
migration  between  carbonate  ions  in  the  cluster  of  (Ei2C03)8.  The  trajectory  analysis  and  the 
activation  energy  strongly  support  that  the  proton  transfer  prefers  the  inter-carbonate  pathway 
and  suggest  that  static  DET  methods  with  a  limited  size  cluster  model  in  the  gas  phase  has 
uncertainties  in  describing  the  state  of  molten  carbonate  and  that  further  refinement  on  calculated 
energies  with  extracted  local  structure  is  necessary.  Such  treatment  can  also  be  applied  in 
treating  similar  molecular  systems  in  future. 


Fig.  4:  MSD  for  proton  transfer  in  the  cluster  of  (LiiCOsjs  simulated  at  1100  K,  1200  K, 
and  1300  K  for  3  ps. 
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Fig.  5:  Arrhenius  plot  for  proton  transfer  in  the  cluster  of  (Li2C03)8 
A3.3:  Mechanism  of  proton  transfer  in  MC 

Due  to  the  faet  that  the  meehanism  of  proton  migration  in  MC  is  eomplieated,  both  statie  DFT 
and  FPMD  caleulations  at  finite  temperatures  have  been  performed  in  the  current  study.  The 
results  from  DFT  calculations  in  the  gas  phase  have  provided  detailed  information  on  structural 
evolution,  charge  transfer  and  energetics,  while  trajectories,  MSD  at  different  finite  temperatures 
and  activation  energy  by  Arrhenius  plot  were  obtained  from  FPMD.  All  theoretical  evidences 
from  multiple  perspectives  have  supported  that  the  inter-carbonate  pathway  is  dominant  for  the 
proton  migration  in  molten  carbonate.  However,  the  intra-carbonate  pathway  is  still  possible  to 
serve  as  an  intermediate  step  in  the  whole  process.  In  addition,  the  combination  of  the  static  and 
dynamic  DFT  methods  in  this  study  has  shown  great  advantages  in  treating  such  ionic  liquid 
systems. 

It  is  interesting  to  note  that  the  degree  of  geometry  change  of  the  carbonate  ions  at  TS  is  slightly 
different  from  the  intra  to  inter-carbonate  pathway.  In  Fig.  1(a),  the  changes  of  the  03  charge 
and  03-C  bond  distance  are  significantly  smaller  than  those  for  other  atoms  and  bonds  in  the 
carbonate  ion,  indicating  03  is  not  actively  participated  in  the  proton  migration  process. 
However,  both  01-C  and  02-C  bonds  change  significantly  during  the  process.  One  bond 
increases  and  the  other  will  decrease  accordingly.  The  Ol-C-02  bond  is  bent  to  106.6°  to 
accommodate  the  proton  migration.  This  is  consistent  with  the  Lewis  diagram  in  Fig.  6(a).  With 
respect  to  the  inter-carbonate  path,  all  atoms  in  both  carbonate  ions  have  experienced  some  small 
changes  in  atomic  charge  and  bond  distance,  showing  they  are  all  actively  involved  in  the 
process.  The  O-C-0  bond  angles  remain  almost  unchanged  (-120°)  in  the  whole  process.  The 
bond  forming  and  breaking  process  can  be  described  by  the  Lewis  diagram  in  Fig.  6(b). 
Theoretically,  if  a  reaction  experiences  more  electron  delocalization  and  less  geometry  distortion 
at  TS,  it  should  have  a  lower  energy  barrier.  From  this  chemical  aspect,  the  inter-carbonate  path 
should  be  more  favored  by  reaction  kinetics. 

For  the  energy  barriers  from  static  DFT  calculations,  they  are  obviously  overestimated.  This 
uncertainty  is  probably  due  to  the  structural  variation  of  the  ionic  cluster,  especially  the  part  not 
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involved  in  the  proton  migration  in  the  optimization  process.  The  value  of  energy  barrier  itself  is 
small,  but  the  total  energy  of  the  cluster  is  large.  To  correct  the  errors,  we  have  extracted  the 
small  local  structure  out  of  the  whole  cluster,  which  allows  us  to  catch  the  energy  change  due  to 
the  proton  migration  and  ignores  the  effect  on  the  energy  from  the  ions  not  involved  in  the 
process.  The  energy  barrier  was  reduced  from  46.8  to  17.5  kcal/mol  for  the  intra-carbonate  path 
and  from  20.8  to  8.0  kcal/mol  for  the  inter-carbonate  one,  respectively.  Such  large  changes  imply 
the  challenges  in  treating  such  an  ionic  liquid  of  molten  salt  state  using  gas  phase  DFT  and 
cluster  models.  The  correction  method  used  here  has  been  very  effective  and  is  easy  to  apply, 
and  further  verifications  will  be  carried  out  in  the  future.  The  final  energy  barrier  from  DFT  is 
8.0  kcal/mol,  which  is  in  a  good  agreement  with  the  activation  energy  of  7.5  kcal/mol  from 
FPMD  simulations.  The  activation  energy  was  measured  to  be  7.6  kcal/mol  using  a  BZY/MC 
composition  electrolyte  at  600  °C  and  MC  was  considered  as  the  main  conducting  media  for 
protons  Based  on  those  values,  we  can  conclude  that  the  energy  barrier  for  the  migration  of 
proton  in  molten  lithium  carbonate  salt  is  in  the  range  of  7. 5-8.0  kcal/mol  and  excellent 
consistency  between  DFT  modeling  and  experimental  measurement  was  observed. 
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Fig.  6:  Lewis  Diagram  for  the  Intra-  (a)  and  Inter- Carbonate  (b)  Migration  of  Proton. 


A4:  Conclusions  for  Part  A 

The  pathway  and  energy  barrier  of  proton  migration  in  the  lithium  molten  carbonate  are 
predicted  by  static  and  dynamic  density  functional  theory  calculations.  The  results  show  that  the 
proton  migration  prefers  the  inter-carbonate  to  the  intra-carbonate  pathway.  The  migration 
energy  barrier  of  8.0  kcal/mol  for  local  structure  extracted  from  the  (Li2C03)8  cluster  is 
consistent  with  the  7.8  kcal/mol  obtained  from  ab  initio  molecular  dynamic  simulations.  Also, 
this  result  is  excellent  in  line  with  the  experimental  value  of  7.6  kcal/mol  observed  in  a  MC/BZY 
composition  at  600  °C.  The  analyses  of  atomic  charges,  bond  distances,  bond  angles,  vibrational 
frequencies,  and  molecular  orbitals  all  indicate  that  the  proton  migration  in  the  lithium  molten 
carbonate  is  a  synergetic  process.  The  trajectories  of  proton,  carbon  and  oxygen  suggest  that  the 
oxygen  atoms  can  adjust  their  positions  during  the  proton  migration  between  the  carbonate  ions 
in  the  cluster  of  (Li2C03)8.  The  MSDs  of  proton  simulations  at  temperatures  of  1100  K,  1200  K, 
and  1300  K  for  3  ps  feature  a  linear  increase  with  time,  implying  a  fast  proton  migration  in  the 
lithium  molten  carbonate.  Moreover,  the  mechanism  of  proton  migration  can  be  explained  by  the 
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Lewis  electron  dot  diagram.  Overall,  the  calculated  results  have  an  excellent  agreement  with  our 
previous  experimental  results  and  indicate  that  molten  carbonate  can  provide  a  low  energy 
barrier  channel  for  proton  conduction  in  IT-SOFCs.  The  kinetics  and  mechanism  of  such  process 
are  well  elucidated  by  DFT  calculations.  In  addition,  our  calculations  indicate  that  the 
combination  of  the  static  and  dynamic  DFT  methods  is  of  great  advantages  in  treating  such  ionic 
liquid  systems  and  can  give  reliable  results. 
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PART  B:  NEW  MC/BZY  HYBRID  PROTON  CONDUCTOR 

Bl:  Experimental  procedures 

Bl.l:  Fabrication  of  MC/BZY  Electrolyte 

The  BZY  (BaZro.sYo.iOs-s)  electrolyte  was  synthesized  via  conventional  solid-state  reaction. 
The  staring  materials  were  BaCOa  (99.8%,  Alfa  Aesar),  ZrOa  (99.7%,  Alfa  Aesar)  and  Y2O3 
(99.9%,  Alfa  Aesar).  Stoichiometric  amount  of  the  starting  materials  were  first  weighed, 
intimately  mixed  and  pelletized,  followed  by  firing  at  1100°C  for  5  h.  The  calcined  pellets  were 
then  broken  into  powders  and  calcined  again  at  1300°C  for  another  5  h.  After  ball  milling,  the 
powders  were  then  blended  with  ~10  wt%  of  carbon -black  as  a  pore  former,  pelletized  and 
finally  sintered  at  ISOO^C  for  5  h.  Thus  made  samples  were  porous  and  ready  for  molten 
carbonate  infiltration. 

Two  molten  carbonate  systems  were  selected  for  this  study:  Li2C03-38mol%  K2CO3 
(denoted  as  LK)  and  Li2C03-  48mol%Na2C03  (denoted  as  LN);  these  compositions  reflect  a 
eutectics.  They  were  first  melted  and  homogenized  in  an  AI2O3  crucible  at  650°C  for  2  h.  The 
porous  BZY  pellet  pre-loaded  in  a  silver  basket  were  then  slowly  immersed  into  the  MC.  After  a 
roughly  2-h  soaking,  they  were  gradually  pulled  out  of  the  melt  and  suspended  above  the 
crucible  before  slowly  cooling  the  finance  to  room  temperature.  Thus  fabricated  MC-BZY 
composites  were  slightly  polished  to  remove  the  residual  carbonate  on  the  surface.  The  samples 
were  then  ready  for  microstructural  and  electrical  characterizations. 

B1.2:  Characterization  of  MC/BZY 

The  phase  purity  of  the  BZY  powder  and  MC-BZY  composites  were  examined  by  X-ray 
diffraction  (XRD)  using  (D/max-A,  Rigaku,  Japan),  while  microstructures  of  the  BZY  as  well 
MC-BZY  composites  were  captured  by  a  field  emission  scanning  electron  microscopy  (FE-SEM, 
Zeiss  Ultra). 

The  ionic  conductivity  of  the  MC-BZY  composites  was  assessed  by  electrochemical 
impedance  spectroscopy  (EIS)  on  a  symmetrical  cell  configuration  consisting  of  the  MC-BZY 
composite  sandwiched  by  two  identical  silver-based  electrodes  and  current  collectors.  The  EIS 
measurements  were  carried  out  using  an  electrochemical  workstation  (IM6,  Zahner)  under  open- 
circuit  conditions  within  a  frequency  range  of  10^  to  0.01  Hz  and  under  a  10  mV  of  stimulus  AC 
amplitude.  The  ionic  resistance  was  extracted  from  the  high-frequency  intersection  of  the 
spectrum  with  the  real-axis.  The  temperature  was  varied  from  400  to  650°C,  covering  the  solid- 
to-liquid  transition  of  the  carbonate  phase  (the  melting  temperature  of  the  carbonate  phase  is 
~490‘’C).  To  understand  the  nature  of  ionic  conduction,  the  ionic  conductivity  was  also  measured 
in  both  oxidizing  (air)  and  reducing  atmospheres  (H2-N2  mixture)  with  a  range  of  varied  partial 
pressures  of  H2O  and  H2.  The  different  partial  pressures  of  H2O  were  generated  by  a  homemade 
water  bubbler  system  and  monitored  by  an  on-line  humidity  sensor  (Vaisala  model  332).  To 
verify  the  involvement  of  protons  in  ionic  conduction,  D2O  was  also  used  as  an  alternative  to 
H2O.  In  addition,  the  effects  of  MC  loading  and  MC-type  on  ionic  conductivity  were  also  studied. 

B1.3:  DFT  modeling  of  proton  transfer  in  molten  carbonates 
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The  energetics  of  proton  transfer  in  a  MC  was  analyzed  by  Gaussian09  suite  of  quantum 
chemistry  programs  [10].  The  disordered  MC  phase  was  simulated  by  a  cluster  of  (Li2C03)8.  The 
B3LYP  hybrid  functional  [11,  12]  in  combination  with  the  basis  set  of  3-21G(d,  p)  were  used  to 
calculate  the  energetics  of  IT^-transfer  in  [(Li2C03)8H]'^  cluster. 


B2:  Results  and  Discussion 

B2.1:  XRD  patterns  and  microstructures 

The  XRD  patterns  of  BZY  synthesized  at  1500°C  for  5  h  is  shown  in  Fig.l  (a),  indicating  two 
cubic  phases,  which  have  been  previously  reported  [13].  The  reactivity  between  BZY  and  LN 
was  also  studied  by  firing  a  mixture  of  BZY  and  LN  (1:1  vol%)  at  650°C  for  2h.  The  XRD 
patterns  are  shown  in  Fig.l  (b),  indicating  no  phases  formed  other  than  the  original  BZY  and  LN, 
which  is  consistent  with  that  reported  in  Ref.  [9].  Rietveld  refinements  of  XRD  patterns  resulting 
from  the  pure  BZY  and  BZY-LN  indicated  that  the  mole  ratio  of  the  two  cubic  BaZr03  phases 
was  consistently  59:41.  The  two  BaZr03  cubic  phases  exhibited  a  lattice  parameter  of 

o  o 

a=4.2021A  and  a=4.2366A,  respectively.  For  the  composite  sample,  the  lattice  parameters  for  the 
two  BaZr03  cubic  phases  were  a=4.1999A  and  a=4.2358A,  further  suggesting  no  reaction 
occurred  between  BZY  and  molten  carbonate. 

The  microstructures  of  BZY  and  MC-BZY  composite  are  shown  in  Fig.2  (a)  and  (b), 
respectively.  It  is  evident  that  BZY  even  after  being  sintered  at  1500°C  contains  poorly 
connected  grains,  indicating  the  refractory  nature  of  BZY  materials.  By  contrast,  Fig.2  (c)  and  (d) 
of  the  MC-BZY  composite  shows  a  dense  microstructure  with  the  MC  phase  filling  the  pores. 


Fig.l  XRD  pattern  of  (a)  BZY  sintered  at  1500“C  for  5h  and  (b)  MC-BZY  fired  at  650“C. 


16 


PROPOSAL  NUMBER:  62940-CH-REP,  FINAL  REPORT  2013-2017 


Fig.2  SEM  images  of  (a)  surface  and  (b)  cross-section  of  BZY  sintered  at  1500"C;  (c) 
surface  and  (d)  cross-section  of  a  LN-BZY  composite  prepared  at  650“C. 

B2.2:  Conductivity  vs  MC  loading 

The  effect  of  MC  loading  on  the  ionic  conductivity  of  a  LK-BZY  is  shown  in  Fig. 3.  Since  the 
MC  loading  is  determined  by  the  initial  porosity  in  BZY,  different  contents  of  carbon-black  pore 
former  were  used  during  preparation  of  porous  BZY.  Fig.3  clearly  indicates  that  higher  MC 
loading  (or  higher  initial  porosity  in  BZY)  results  in  a  higher  ionic  conductivity  over  the  entire 
temperature  range  tested.  This  is  an  expected  result  as  the  primary  conductor  in  MC-BZY  is  the 
MC  phase,  which  is  also  evidenced  by  the  activation  energy  Ea=0.29-0.30  eV  close  to  that  of 
pure  molten  carbonate.  In  practice,  however,  the  high  conductivity  of  the  MC  phase  could  be 
limited  by  the  amount  of  MC  that  can  be  retained  within  the  BZY  skeleton,  analogous  to  the  case 
of  MC  in  LiA102  matrix  used  in  commercial  molten  carbonate  fuel  cells. 
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Fig.3  Ionic  conductivity  vs  temperature  in 
LK-BZY  composites  at  different  MC 
loadings 
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B2.3:  Conductivity  vs  MC  type 

The  effect  of  MC  type,  i.e.  LK  vj'  LN,  on  ionic  conductivity  of  MC-BZY  is  shown  in  Fig.4.  To 
ensure  a  fair  comparison,  the  MC  loading  of  the  two  samples  was  fixed  at  48  vol%.  Clearly,  LN- 
BZY  has  a  higher  conductivity  than  LK-BZY  at  t>500°C.  This  observation  is  consistent  with  the 
literature  showing  higher  conductivity  for  LN  than  LK  system  [14].  In  addition,  Ea=0.30  and 
0.40  eV  at  above  the  melting  temperature  for  the  LN-BZY  and  LK-BZY  system,  respectively, 
are  comparable  to  the  literature  values  [14].  The  LN  system  was  selected  for  the  subsequent 
electrical  characterization  because  of  its  high  conductivity. 


Fig.4  Conductivity  comparison  of  a  LK- 
BZY  and  LN-BZY  vs  temperature  in 
3%H20-air. 
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B2.4:  Conductivity  vs  atmosphere 

Fig. 5  shows  the  conductivity  of  a  LN-BZY  with  a  MC  loading  of  57  vol%  measured  in  both 
oxidizing  and  reducing  atmospheres.  It  is  evident  that  conductivity  measured  in  reducing 
3%H20-(5%H2-N2)  is  higher  than  that  measured  in  oxidizing  3%H20-air  over  the  entire 
temperature  range.  The  result  suggests  that  the  enhanced  conductivity  by  H2  be  associated  with 
the  MC  phase  since  the  BZY’s  conductivity  above  the  melting  temperature  is  negligible 
compared  to  the  MC  phase.  The  almost  identical  activation  energy  between  the  LN-BZY  and 
pure  LN  shown  in  Fig. 5  supports  this  assertion. 
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Fig.5  Conductivity  vs  temperature  of  a 
LN-BZY  composite  measured  in  3%H20- 
air  and  3%H20-(5%H2-N2) 
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B2.5:  Conductivity  vs  Ph^o 

To  further  investigate  the  feasibility  of  proton  conduction,  the  conductivity  of  a  LN-BZY 
composite  with  a  MC  loading  of  57  vol%  was  measured  as  a  function  of  partial  pressure  of  H2O 
( PH20)  at  bOO^C;  the  results  are  shown  in  Fig.6.  First,  the  conductivity  in  reducing  5%H2-N2  is 
nearly  twice  as  high  as  that  in  air  under  each  PH20,  which  agrees  with  Fig. 5  measured 
under  Pf^2O=0-03  atm.  Second,  the  conductivity  increases  withPf^^o-  Third,  the  dependence  of 
conductivity  on  is  almost  the  same  for  the  two  gas  compositions.  These  observations  again 
suggest  that  the  proton  conduction  be  associated  with  the  MC  phase  and  protons  be  originated 
from  either  H2O  or  H2  or  both. 

From  a  defect  chemistry  perspective,  the  enhanced  proton  conductivity  by  H2  and  H2O  can  be 
understood  by  the  following  reaction  (in  Kroger-Vink  notation): 

H20(g)  +  Vo  +  (C03)^03  =  (OH)o  +  (HC03)oo3  (D 

Here  Vq*  denotes  oxygen  vacancies  in  BZY,  serving  as  the  charge  compensation  carriers;  (OHq) 
represents  hydroxyl  defect  residing  on  oxygen  lattice  of  BZY;  (C03)co3  is  a  regular  carbonate-ion; 
(HC03)co3  is  proton  conducting  species  residing  on  carbonate-ion.  From  reaction  (1),  it 

is  evident  that  increasing  partial  pressure  of  H2O  can  increase  the  concentrations  of  (OHo)and 
(HC03)co3  ,  ultimately  leading  to  an  increase  in  conductivity  by  protons.  It  is  also  inferred  from 
Fig.6  that  (HC03)co3  residing  on  carbonate-ion  is  the  dominant  charge  carrier  for  proton 
conduction.  It  is  worth  mentioning  that  reaction  (I)  implicates  a  synergetic  proton  transfer 
between  BZY  and  MC,  in  which  bothY”  and(C03)co3Play  an  important  role.  The  details  on  the 
proton  transfer  mechanism  will  be  discussed  in  the  following. 


Fig.6  Conductivity  of  a  LN-BZY 
composite  vs  Puzo  at  600"C. 


B2.6:  Conductivity  vs  Ph^ 

To  understand  the  enhanced  conductivity  by  H2,  the  conductivity  of  LN-BZY  composite  was 
measured  over  a  range  of  partial  pressure  of  H2  (Ph^)  at  600°C;  the  result  is  shown  in  Fig.7.  Note 
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that  the  conductivity  shown  in  Fig.7  is  lower  than  that  in  Fig. 6.  This  is  caused  by  a  lower  MC 
loading  in  the  sample,  viz.  48  vol%  for  Fig.7  vs  57  vol%  for  Fig. 6.  The  defect  chemistry  reaction 
leading  to  the  increased  concentrations  of  (OHq)  andlHCOjl^o, :  and  thus  the  conductivity  can 
be  expressed  by: 

H,  +  y"  +  +  (CO.Tco,  =  (OHYo  +  (HCO^Yco^  (2) 

here  Og  is  the  regular  oxygen  lattice  in  BZY.  Again,  ( HCO^  is  the  dominant  charge  carrier 

over  because  MC  is  the  primary  ionic  conductor  at  600°C.  The  synergy  of  proton  transfer 

between  BZY  and  MC  phase  is  further  discussed  in  the  following. 


0.2  0.4  0.6  0.8 

Partial  pressure  of  H^,  atm 


Fig.7  Conductivity  of  a  LN-BZY 
composite  vs  Ph2  at  600‘'C 


B2.7:  Conductivity  in  D2O 


To  verify  whether  H  in  H2O  is  one  of  the  proton  sources,  we  measured  conductivity  in  the 
presence  of  D2O.  Fig. 8  compares  the  conductivities  measured  in  air  saturated  with  3%H20  and 
3%D20.  The  conductivity  difference  for  the  two  gases  is  small  but  visible,  suggesting  the  nature 
of  proton  conduction.  As  aforementioned,  proton  species  (HCO^Ygo^  *^he  primary  source  for  the 


observed  enhanced  conductivity  in  both  H2O  and  H2  containing  atmospheres.  The  close  Ea  to 
pure  LN  shown  in  Fig. 8  further  indicates  the  dominance  of  ionic  conduction  by  the  molten 
carbonate  phase. 


Fig.8  Conductivity  of  a  LN-BZY  vs 
temperature  in  3%H20-air  and  in 
3%D20-air. 


20 


PROPOSAL  NUMBER:  62940-CH-REP,  FINAL  REPORT  2013-2017 


B2.8:  DFT  modeling  of  proton  transfer  in  molten  carbonates 

To  understand  the  theoretical  feasibility  of  proton  conduction  in  molten  carbonates,  we 
performed  a  first-principles  DFT  calculation  of  the  energetics  of  proton  transfer  in  a  molten 
carbonate  phase.  Fig.9  shows  the  energetics  of  proton  transfer  intrinsic  reaction  coordinate  (IRC) 
across  a  cluster  of  LiiCOa  (Fig.9  (a)),  a  surrogate  for  molten  carbonate  for  the  DFT  calculations. 
The  marked  in  green  color  bonds  with  05  as  the  reactant  and  bonds  with  020  as  the  product. 

The  IRC  calculations  specify  that  the  first  rotates  counter-clockwise  along  the  05  to  reach  the 

0 

TS.  During  this  process,  the  bond  between  H  and  05  is  slightly  shortened  from  1.070  A  to 

0 

1.030  A  with  an  energy  barrier  of  8.9  kcal/mol.  The  conversion  from  TS  into  product  is  rather 
fast.  The  separates  from  05  and  moves  towards  020.  Meanwhile,  the  bond  length  of  H'^-05 

o  o  o 

increases  from  1.030  A  to  1.267  A,  whereas  the  bond  length  of  H  -020  reduces  from  2.933  A  to 

0 

1.174  A,  thus  completing  a  single  H  -transfer  process.  In  addition,  it  is  also  seen  from  Fig.9  (b) 
that  53.1  kJ/mole  energy  is  needed  to  overcome  for  a  reversible  transition  from  product  to 
reactant,  which  further  implies  that  continuous  IT^-transfer  in  (Li2C03)8  cluster  or  molten 
carbonate  is  feasible. 


Fig.9  The  energetics  of  proton  transfer  in  cluster  of  (Li2C03)8;  (a)  cluster  structure;  (b)  the 

relative  energies  vs  IRC.  Unit:  kcal/mole 
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B2.9:  A  synergetic  proton  transfer  mechanism 

As  indicated  by  reaction  (1)  and  (2),  the  creation  of  proton  species  replies  upon  the 

presence  of  V"  and  Og  in  the  BZY.  In  other  words,  proton  transfer  across  the  composite  is 
synergetic  between  BZY  and  MC.  To  illustrate  the  mechanism,  Fig.  10  shows  a  schematic  of  the 
synergetic  proton  transfer  in  the  atmospheres  containing  H2O  or  H2  or  both.  It  is  worth 
reemphasizing  that  (HCO^Ygg^  is  considered  the  primary  proton  conducting  species 
accompanying  the  major  carbonate-ion  conduction  in  the  MC-BZY  composite  electrolyte. 


H2O  H2 


Fig.lO  Schematic  of  synergetic  proton  conduction  in  MC-BZY  composite  electrolyte 
B3:  Conclusion  and  Summary  for  Part  B 

In  summary,  electrical  properties  of  the  MC-BZY  composite  electrolyte  was  methodologically 
evaluated  to  understand  the  enhanced  conductivity  and  improved  SOFC  performance.  It  was 
shown  that  the  composite  can  be  made  into  a  dense  membrane  at  bSO^C  by  filling  the  pores  in 
BZY  with  MC  phase.  The  ionic  conductivity  of  the  composite  above  500°C  was  affected  by  the 
loading  and  type  of  MC.  More  important,  the  sample  exposed  in  H2-containing  atmospheres 
exhibited  nearly  twice  the  conductivity  of  the  one  in  air.  The  dependences  of  conductivity  on 
partial  pressures  of  H2  and  H2O,  as  well  as  the  type  of  H  (H  vs  D)  suggested  that  proton 
conduction  take  places,  primarily  in  the  MC  phase.  A  defect  chemistry  model  was  presented  to 
explain  the  synergetic  proton  transfer  between  BZY  and  MC  phases.  A  preliminary  first- 
principles  DFT  calculation  further  confirmed  that  proton  transfer  in  molten  carbonate  was  indeed 
a  low-energy  barrier  process. 
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PART  C:  PROTON  MIGRATION  AT  THE  BZY  AND  MC  INTERFACE 
Cl:  Introduction 

A  focus  of  the  recent  development  of  solid  oxide  fuel  cell  (SOFC)  technology  towards 
commercialization  is  to  lower  SOFC’s  operating  temperature  from  the  current  700-800°C  to 
below  bOO^C.  One  key  to  realize  such  reduced  temperature  SOFCs  is  to  discover  and  employ  a 
high-conductivity  solid-oxide  electrolyte.  One  promising  class  of  solid  oxide  electrolytes  is 
BaMOa  (M=Zr,  Ce)  based  perovskite  protons  (H"^)  conductors  Unlike  conventional  intrinsic 
H'^-conductors  such  as  H3PO4,  these  oxide -based  H'^-conductors  strongly  rely  on  the  oxygen- 
vacancy  ( V” )  and  its  interactions  with  H2O  to  produce  IF^-carrying  species  of  OHq  ,  which  is 
known  as  the  “Grotthuss  mechanism”  Since  the  hydration  reaction  is  exothermic,  lower 
temperatures  favor  the  formation  of  OHq  ,  thus  proton  conductivity.  This  fundamental  fact  is 
interestingly  synergetic  with  low-temperature  SOFCs,  making  BaMOa-based  materials  a 
promising  low -temperature  solid  electrolyte  candidate. 

BaZro.8Yo.2O2.95  (BZY)  is  a  representative  of  well-studied  oxide-based  proton  conductors  with 
reasonably  high  conductivity  Compared  with  its  counterpart  BaCeOa,  BZY  attracts  much 
more  attention  due  to  its  chemical  stability  in  C02-containing  atmospheres.  However,  it  is  hard 
to  sinter  it  into  a  dense  microstructure,  rendering  difficulties  for  practical  use.  Recently,  we 
showed  that  molten  carbonate  (MC)  salts  infiltrated  into  BZY  can  help  produce  a  dense 
microstructure  More  interestingly,  such  a  dual -phase  electrolyte  exhibited  enhanced  proton 
conductivity.  Similar  results  have  also  been  observed  by  other  groups  with  excellent  SOFC 
performance  In  an  effort  to  understand  this  phenomenon,  we  proposed  that  OHq  derived 
fromV”  transfers  to  the  neighboring  carbonate-ion  (COjI^q^  ,  forming  a  new  proton  species 
(HCOjIqo^  MC  phase  **’'^*,  through  the  following  defect  reaction  (in  Kroger-Vink  notation): 

H20(g)  +  Vo  +  (C03)^oo3  =  (OH)o  +  (HC03)*  03  (1) 

A  subsequent  DFT  modeling  revealed  that  proton  migration  in  the  form  of  (HCOjI^q  within  a 

MC  phase  is  a  very  facile  process  with  low  energy  barrier  However,  how  H"^  is  transferred 
across  the  interface  of  BZY/MC  as  suggested  by  reaction  (1)  remains  ambiguous  at  this  point 
despite  the  fact  that  a  range  of  computational  modeling  work  on  H"^  migration  and  stability  in 
BaMOa-based  proton  conductors  has  been  reported  in  the  open  literature 

In  the  present  work,  we  aim  to  understand  the  energetics  of  proton  migration  across  the  BZY/MC 
interface  from  a  DFT  perspective.  Pathways  and  energetics  of  proton  migration  at  the  Zr02- 
terminated  surface  (100)  of  a  pure  BZ  (BZ:  BaZrOa)  and  a  BZ  with  MC  are  particularly 
computed  and  compared.  The  results  clearly  suggest  that  the  presence  of  MC  can  significantly 
enhance  the  H"^  migration  across  the  surface  of  BZ  (and  BZY). 

C2:  Computational  method  and  model 

All  DFT  calculations  conducted  in  this  study  were  performed  using  a  Vienna  ab  initio  simulation 
package  (VASP)  with  the  projector  augmented  wave  (PAW)  approach  The  exchange- 
correlation  term  of  electrons  was  described  by  the  Perdew-Burke-Ernzerhof  (PBE)  functional 
Eor  all  surface  calculations,  a  4x4x1  Monkhorst-Pack  sampling  of  the  Brillouin  zone  was  used 
A  vacuum  region  of  15  A  in  the  z-axis  direction  was  created  to  prevent  any  interaction  between 
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cells.  The  O  2s^2p'^,  C  2s^2p^,  Zr  4s^4p^5s^4d^,  4s^4p^5s^4d'  of  Y,  and  Ba  5s^5p^6s^  were  taken 
as  the  valence  electron  configurations  in  the  calculations.  To  ensure  a  high  accuracy,  a  cutoff 
energy  of  500  eV  was  chosen  for  the  plane  waves.  During  the  structural  optimization,  all  atoms 
were  fully  relaxed  with  a  force  less  than  0.02  eV/A.  Furthermore,  the  D2  method  of  Grimme 
was  applied  to  calculate  the  dispersion  interaction  between  carbonate-ion  and  surface,  in  which 

the  cutoff  radius  for  pair  interactions,  global  scaling  factor  optimized  at  the  PBE  and  damping 

0 

parameter  were  taken  as  15.0  A,  0.75  and  20.0,  respectively.  The  saddle  points  and  migration 
pathways  with  the  minimum  energy  were  sought  by  the  nudged  elastic  band  (NEB)  method 
Eor  all  the  calculations,  spin-polarized  and  dipole  corrections  have  also  been  taken  into  account. 

Calculations  on  energetics  and  pathways  of  proton  migration  were  conducted  on  an  8 -layer 
asymmetrical  slab  at  the  Zr02-terminated  (100)  surface.  The  (100)  surface  was  selected  because 
it  has  been  previously  determined  to  be  the  most  stable  surface  for  BaZrOs  An 

asymmetrical  slab  was  chosen  to  maintain  a  stoichiometric  and  neutral  unit  cell,  which  is  an 
acceptable  approach  in  the  study  of  similar  systems.  Eor  example,  Evarestov  et  al  showed  a 
negligible  difference  in  surface  energy  between  an  asymmetric  and  symmetric  EaMnOs  (110) 
surface  slab.  Eiu  et  al  used  an  8-layer  asymmetric  EaMnO -terminated  EaMnOa  (110)  surface 
slab  in  Ref  30.  Eurthermore,  previous  studies  have  shown  that  a  total  of  8  layers  are  sufficient  to 
simulate  the  surface  processes  on  BaZrOs  In  the  current  model,  the  bottom  4  layers 

were  fixed  to  represent  the  bulk  of  BaZrOa,  whereas  the  top  4  layers  were  relaxed  and  considered 
as  the  domain  for  adsorbed  and  COa^'.  On  the  other  hand,  with  a  ZrOa-terminated  surface, 
can  be  modeled  to  migrate  along  the  edge  of  a  ZrOe  octahedron,  which  has  been  previously 
determined  as  a  favorable  pathway 


C3:  Results  and  Discussion 

C3.1:  Adsorption  of  proton  on  the  Zr02-terniinated  surface 

We  started  calculations  by  first  optimizing  the  bulk  structure  with  a  6x6x6  k-point  mesh  using 
the  computational  methods  described  above.  The  calculated  lattice  constant  of  a  cubic  BaZrOa  is 
4.25  A,  which  is  in  good  agreement  with  the  experimental  4.20  A  and  theoretical  4.24  A  . 
The  adsorption  behavior  of  a  single  proton  on  the  ZrOa-terminated  surface  was  then  studied.  A 
total  of  five  adsorption  sites  were  considered  and  are  schematically  shown  in  Eig.l:  (a)  on  the  top 
of  O  atom  (upright).  To,  (b)  on  the  top  of  O  atom  (tilted  to  Zr  atom),  To-zr,  (c)  on  the  top  of  O 
atom  (tilted  to  hollow  site),  To-hoiiow,  (d)  on  the  top  of  Zr  atom,  Tzr,  and  (e)  on  the  top  of  hollow 
site,  Thoiiow.  After  a  full  relaxation,  the  To-hoiiow  configuration  (c)  is  calculated  to  be  the  most 
favorable  adsorption  site  with  the  lowest  energy,  indicating  that  the  proton  prefers  to  be  on  the 
top  of  O  and  tilted  to  the  hollow  site.  This  finding  is  in  agreement  with  previous  understanding 
on  the  formation  of  proton  defects,  where  water  from  the  gas  phase  dissociates  into  a  hydroxide 
ion  and  a  proton.  The  former  will  then  fill  an  oxygen  vacancy,  while  the  latter  forms  a  covalent 
bond  with  the  lattice  oxygen  Moreover,  Eig.l  (c)  also  indicates  that  the  lattice  oxygen  deviates 
from  its  lattice  position  to  the  hollow  site  and  that  the  proton  tilts  to  the  hollow  site,  shortening 
the  distance  between  proton  and  its  neighboring  oxygen.  Therefore,  this  geometry  is  in  favor  of 
the  proton  migration  from  an  O  to  its  adjacent  O,  so  will  reduce  the  energy  barrier  of  proton 
transfer. 
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(a)  Etot  =  -660.423  eVv  (b)  E,ot  =  -659.050  eV^ 


(d)  E,ot  =  -659.326  eVv  (e)  Etot  =  -659.287  eV^ 


(c)  Etot  =  -660.746  eV^ 

• 

Ba  Zr  O 
X.M - 

▼ 


Fig.l  Optimized  structures  (top  view)  and  total  energy  of  proton  adsorbed  on  the  ZrOi- 
terminated  surface,  (a)  on  the  top  of  O  atom  (upright),  To,  (b)  on  the  top  of  O  atom  (tilted 
to  Zr  atom),  To-zr,  (c)  on  the  top  of  O  atom  (tilted  to  hollow  side),  To-hoiiow,  (d)  on  the  top  of 
Zr  atom,  Tzr,  and  (e)  on  the  top  of  hollow  site,  T  hollow. 

C3.2:  Proton  migration  on  the  Zr02-terminated  surface 

Next,  selected  migration  pathways  and  their  corresponding  energy  barriers  for  the  proton 
migration  on  the  Zr02-terminated  (100)  surface  were  studied.  Fig.2  shows  the  migration 
pathways  and  energy  profiles  of  along  X-direction  (a  to  b),  Y-direction  (c  to  d),  and  XY 
direction  (e).  For  the  X-direction,  two  independent  migration  pathways  are  possible:  linear  path  1 
2  (labeled  as  “Path-1”)  with  the  shortest  migration  distance  and  curved  path  from  1  ^  3  ^  4 
2  (labeled  as  Path-2).  These  two  migration  pathways  and  the  corresponding  energetics  are 
shown  in  Fig.2  (a)  and  (b),  respectively.  “Path-1”  includes  one-step  migration  and  experiences  an 
energy  barrier  of  1.660  eV.  “Path-2”  consists  of  three  steps,  marked  as  I,  II,  and  III.  The  step  I 
and  III  are  related  to  transfer  from  an  oxygen  toward  a  neighboring  oxygen  with  the  bond 
breaking  and  reforming  of  O-H"^.  The  step  II  involves  only  rotation  around  an  oxygen  with  the 
O-H"^  bond  connected.  The  energy  barriers  of  transfer  and  rotation  are  0.567  eV  and  0.326  eV, 
respectively.  Therefore,  the  overall  energy  barrier  for  “Path-2”  is  0.567  eV,  which  is  in 
agreement  with  the  activation  energy  of  0.69  eV  for  IT^  transfer  observed  in  BaZrOs  mobility 
activation  enthalpies  of  0.41-0.63  eV  for  protonic  defects  in  cubic  perovskite-type  oxides 
activation  energy  of  0.3-0. 6  eV  and  0.45-0.65  eV  for  IT^  migration  in  BZY  electrolytes. 
Evidently,  the  proton  transfer  is  the  rate-determining  step,  which  agrees  well  with  the 
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experimental  and  other  computational  results  Overall,  the  results  show  that  proton 

migration  preferentially  favors  a  curved  over  a  straight  pathway. 


(a).  Reaction  Path  (A) 


(h)»  Reaction  Path  (A) 
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(d).  Reaction  Path  (A) 


Fig.2  Migration  pathways  and  energy  profiles  for  H+  migration  along  X  direction  (a  to  b), 
Y  direction  (c  to  d),  and  XY  direction  (e).  The  migration  pathways  are  marked  with  arrows. 

Fig.2  (c)  and  (d)  show  two  more  different  migration  pathways,  labeled  as  “Path-3”  and  “Path-4”, 
for  migration  along  the  Y-direction.  These  two  pathways  represent  migration  along  a 
linear  and  zigzag  pattern,  respectively.  Fig.2  (c)  shows  that  “Path-3”  follows  migration  sequence 
1  ^  2  ^  3  ^  1 ,  consisting  of  a  rotation  step  (I)  followed  by  a  transfer  step  (II)  via  the  top  of 
hollow  site.  The  former  needs  to  overcome  an  energy  barrier  of  0.326  eV,  while  a  much  higher 
energy  barrier  of  2.187  eV  for  the  latter.  The  high  energy  profile  suggests  that  IT^  migration  via 
this  pathway  be  energy  prohibitive.  Fig.2  (d)  shows  the  “Path-4”  via  1^2^4^3^1, 
consisting  of  a  rotation  step  (I)  followed  by  two  transfer  steps  (II,  III).  The  highest  energy 
barrier  for  IT^  migration  is  0.569  eV,  indicating  that  the  zigzag  pathway  is  the  favorable  route  for 
IT^  migration. 

From  the  symmetry  viewpoint,  we  also  consider  XY-direction  as  the  “Path-5”  for  migration. 
Fig.2  (e)  shows  “Path-5”  following  a  migration  sequence  ofl^2^3^4^1,  which  consists 
of  two  rotation  steps  (I  and  III)  and  two  transfer  steps  (II  and  IV).  The  highest  energy 
barrier  for  this  type  of  migration  is  0.567  eV,  also  indicating  a  favorable  pathway. 

Based  on  the  DFT  analysis  above,  we  now  can  conclude  that  migration  on  the  ZrOi- 
terminated  (100)  surface  of  BaZrOs  favors  a  curved  pathway,  involving  both  transfer  and 
rotation  of  proton  conduction  mechanism. 
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C3.3:  Proton  migration  on  the  ZrOi-terminated  surface  with  Y-doping 


For  practical  use,  Y-doping  on  Zr-site  of  BaZrOa  is  necessary  to  introduce  oxygen  vacancies,  and 
the  latter  is  the  source  of  proton  species  OHq.  The  model  for  the  ZrOa-terminated  (100)  surface 
with  an  yttrium  on  the  Zr-site  is  shown  in  Fig. 3  (a).  We  assumed  that  an  oxygen  vacancy  will  be 
created  next  to  yttrium  and  this  allows  to  examine  its  effect  on  energy  barrier  of  migration. 


o  c 

(a). 

Ba  Zr  Y  O 


(b). 


Fig.3  Migration  pathways  and  energy  profiles  for  H+-migration  along  Y  direction  at  Y- 
doped  Zr02-termination  surface  (100).  The  migration  pathways  are  marked  with  arrows. 

According  to  the  coordination  environment,  two  types  of  O,  denoted  as  A  and  B,  are  identified  in 
the  model.  The  A-type  O  bonds  with  Y  to  form  a  half  YOe  octahedron  on  the  surface,  while  the 
B-type  O  bonds  with  Zr  to  form  another  half  ZrOe  octahedron  on  the  surface.  To  compare  with 
the  energy  barriers  of  migration  on  a  perfect  surface  of  Zr02,  migration  along  Y-direction 
(“Path-4”)  at  the  defective  surface  was  examined.  The  migration  pathway  via 
1^2^3^4^5  includes  four  steps:  step  I  and  IV  represent  rotation  around  A-and  B- 
type  O,  respectively;  step  II  and  III  represent  transfer  between  two  A-type  O  and  between  A- 
and  B-type  O,  respectively.  Energy  barriers  of  the  four  migration  steps  are  given  in  Fig.3  (b), 
with  0.266,  0.533,  0.523,  and  0.340  eV  for  step  I,  II,  III  and  IV,  respectively.  The  total  energy 
barrier  for  IT^  migration  along  this  pathway  is,  therefore,  0.533  eV,  which  is  not  significantly 
different  from  0.567  eV  of  the  perfect  surface,  implying  that  the  effect  of  yttrium  on  the  energy 
barrier  of  migration  on  the  surface  is  negligible.  Therefore,  the  following  DFT  calculations  of 
MC  effect  on  IT^  migration  were  primarily  conducted  on  a  perfect  BZ  surface  with  ZrOi- 
termination. 

C3.4:  Proton  migration  at  the  Zr02-terminated  surface/MC  interface 

To  understand  the  observed  enhanced  proton  conductivity  of  BZY  by  MC,  we  further  calculated 
the  energy  barrier  of  migration  at  the  interface  consisting  of  a  Zr02-terminated  (100)  surface 
and  MC  using  the  method  described  above.  The  migration  along  the  Y-direction  (“Path-3”) 
was  used  as  an  example  for  illustration  purpose.  Fig.4  shows  the  pathway  considered  and  energy 
barriers  calculated  for  IT^  migration,  where  (a)  and  (b)  are  the  side  and  top  view,  respectively, 
and  (c)  is  the  energetics  calculated  for  each  migration  step.  The  pathway  consists  of  four 
individual  steps:  1^2^3^4^1.  Step  I  ( 1  — >  2 )  is  the  rotation  around  an  O  of  the  Zr02- 
terminated  surface;  step  II  ( 2  ^  3 )  is  the  transfer  from  surface  to  COa^',  step  III  ( 3  ^  4 )  is 
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rotation  around  an  O  of  COa^',  and  finally,  step  IV  ( 4  ^  1 )  is  transfer  from  COa^'back  to 
surface.  It  is  interesting  to  see  that  during  the  whole  process  of  migration,  COa^'  links  position 
“2”  to  “1”  like  a  bridge  lying  above  the  hollow  site.  Fig.4  (c)  displays  that  energy  barriers  for 
step  I,  II,  III  and  IV  are  0.326,  0.281,  0.332,  and  0.061  eV,  respectively.  The  highest  energy 
barrier  is  0.332  eV,  suggesting  that  the  step-III  of  rotation  around  O  in  COa^'  is  the  rate- 
limiting  step,  whereas  IT^  transfer  between  surface  and  COa^'  is  relatively  easy  with  a  low 
energy  barrier.  The  energy  barrier  of  0.332  eV  for  IT^  migration  is  in  excellent  agreement  with 
0.33  eV  obtained  from  the  conductivity  measurement  of  a  BZY-MC  composite  Compared  to 
the  results  shown  in  Fig. 3,  the  presence  of  COa^'  on  the  surface  of  BaZrOa  has  significantly 
reduced  the  energy  barrier  for  IT^  migration  from  0.567  eV  to  0.332  eV.  Overall,  the  calculations 
confirm  that  MC  provides  an  alternative  easy  pathway  for  protons  to  travel  along  the  surface  of 
BZY.  The  observed  enhanced  proton  conductivity  in  BZY/MC  composite  is,  therefore,  the  result 
of  a  facile  proton  exchange  between  BZY  and  MC  phase  along  the  interface. 


- 

(a)  side  viewv 


(b)  top  vieww 


o 


Ba  Zr  O  C  H 


(c)  V  Reaction  Path  (A) 


Fig.4  Migration  pathways  and  energy  profiles  (c)  for  proton  migration  along  Y  direction 
with  the  aid  of  molten  carbonate  (MC),  (a)  is  side  view  and  (b)  is  top  view.  The  migration 
pathways  are  marked  with  arrows. 
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C4:  Conclusion 

In  summary,  we  have  for  the  first  time  studied  proton  migration  at  the  interface  of  BaZr(Y)03 
and  molten  carbonate  by  density  functional  theory  calculations.  Proton  migration  along  the  Zr02- 
terminated  (100)  surface  of  BaZrOs  favors  a  curved  path  with  an  energy  barrier  of  0.567  eV.  In 
the  presence  of  carbonate-ion,  the  energy  barrier  is  significantly  reduced  to  0.332  eV,  which  is  in 
excellent  agreement  with  0.33  eV  of  activation  energy  observed  in  conductivity  measurement. 
The  rate-determining  step  is  the  proton  transfer  with  the  absence  of  MC,  whereas  it  becomes  the 
proton  rotation  around  oxygen  in  carbonate-ion  in  the  presence  of  MC.  In  addition,  the  effect  of 
Y-doping  on  the  energy  barrier  for  proton  migration  on  surface  is  insignificant  (0.567  eV  for  BZ 
surface  vs  0.533  eV  for  BZY  surface).  Overall,  the  enhanced  proton  conductivity  in  BZY-MC 
composite  electrolyte  is  attributed  to  the  facile  interfacial  proton  exchange  between  BZY  and 
MC  phase. 
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PART  D:  MOLTEN  CARBONATE  AS  A  NON-METAL  CATALYST  FOR  CO 
OXIDATION 

DI:  Introduction 

Carbon  monoxide  (CO)  is  not  only  toxic  to  human  and  environment  but  also  a  poison  to  many 
metal  catalysts  in  industrial  processes.  The  easiest  way  to  remove  CO  from  other  gases  is  to  turn 
it  into  carbon  dioxide,  a  less  toxic  and  inert  gas.  Precious  metals  have  been  mostly  promoted  for 
such  purpose  with  promising  performances  observed.  For  example,  Bond  et  al.  have  firstly 
demonstrated  Pd/SnOi  as  an  excellent  CO  oxidizing  catalyst  at  the  temperature  of  120-200  °C  [1]. 
More  recently.  Stark  and  Harris  developed  the  catalyst  systems  of  Pd/Sn02  and  Pt/Sn02, 
showing  every  high  conversion  efficiency  even  at  -27  °C.  They  were  successfully  deployed  to 
CO2  laser  systems  [2-3].  Additionally,  other  metal  catalysts  such  as  Au/Mn02  [4-6],  Au/Fe203 
[7-8],  Pt/AhOs  [9-11],  Ru-Pt  [12-13]  and  Ag/Mn02  [14-15]  have  also  been  extensively 
investigated  for  catalyzing  the  oxidation  of  CO  at  ambient  or  low  temperatures.  Significent 
potential  applications  for  such  systems  include  CO2  laser,  auto  exhaust  gas  treatment,  gas 
purification  and  mine  rescue. 

To  reduce  or  remove  the  usage  of  precious  metals  as  they  are  of  high  cost  and  limited  resource 
availability,  many  efforts  have  been  geared  towards  low-cost  transition  metals  and  their  oxides 
[16-18].  For  example,  Wang  et  al.  prepared  a  copper-based  catalyst  of  7.5Cu/10SDC  and 
achieved  a  conversion  of  97%  at  175  °C  [19].  Kang  et  al.  reported  a  catalyst  of  CoOx/Ce02  with 
a  very  low  reaction  temperature  and  good  water  vapor  resistance  [20].  It  was  proposed  that 
highly  dispersed  Co  at  high  valence  state  palys  an  important  role  in  the  process  of  CO  oxidation. 
More  recently,  nanostructures  like  graphene  [21-22],  nanotube  [23-24],  zeolite  [25],  and 
MOF[26-27]  have  also  been  extensively  examined  for  CO  oxidation.  As  we  know,  the  gas 
mixtures  are  at  very  high  temperatures  in  some  cases  such  as  combustion  engine  exhaust  and 
industry  processes,  which  may  cause  changes  of  the  geometrical  and  electronic  structure,  and 
chemical  reactivity  for  the  metal-based  catalysts.  Finally,  the  catalyst  will  be  degraded. 

In  the  past  decades,  molten  carbonates  (MC)  have  been  used  in  the  fuel  cells  as  electrolyte. 
Recently,  they  were  introduced  into  the  cathode  of  solid  oxide  fuel  cells  (SOFC)  as  a  booster  for 
the  reacton  of  oxygen  reduction.  Qin  et  al.  studied  the  dissociative  adsorption  of  oxygen  in  MC 
using  density  functional  theory  (DFT)  modeling  and  found  that  oxygen  exists  in  MC  as 
perocarbonate  (C04^')  [28].  In  C04^',  an  oxygen  atom  is  carried  by  a  carbonate  ion  and 
considered  as  a  strong  oxidizing  agent.  To  prove  the  concept  that  MC  could  be  a  possible  non- 
metal  catalyst  for  oxidization,  we  have  performed  a  combined  experimental  and  computational 
study  of  CO  oxidation  catalyzed  by  MC  in  the  current  report.  To  our  knowledge,  this  is  the  first 
time  to  use  liquid  phase  molten  salt  to  catalyze  the  CO  oxidation  and  will  open  a  door  for 
developing  non-metal  eco-friendly  catalysts  for  different  industrial  and  environmental  processes 
in  the  future. 

D2:  Experimental  and  Computational  Methods 
D2.I:  Catalyst  preparation 

Ternary  eutectic  molten  carbonate  of  Li/Na/K  was  supported  on  AI2O3  (1/8”  pellets,  high  surface 
area,  Alfa  Aesar)  and  used  as  CO  oxidation  catalyst.  Li2C03  (Alfa  Aesar,  99.8%),  Na2C03  (Alfa 
Aesar,  99.9%)  and  K2CO3  (Alfa  Aesar,  99.8%),  the  precursors  of  the  ternary  molten  carbonate, 
were  dissolved  in  deionized  water  in  a  molar  ratio  of  42.5  :  32.5  :  25.0  to  form  a  saturated 
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aqueous  solution.  The  catalyst  (MC/AI2O3)  was  prepared  by  impregnation  of  porous  AI2O3  with 
the  ternary  carbonate  aqueous  solution.  More  specifically,  the  AI2O3  porous  pellets  were  soaked 
in  carbonate  solution  for  approximately  4  hours  in  a  vacuum  chamber  (lOKPa),  then  dried  in  an 
oven  at  80  °C  for  12  hours  and  followed  with  calcinating  at  400  °C  to  remove  the  moisture.  The 
MC  load  was  approximately  1.0%  by  weight. 

D2.2:  Catalyst  characterization 

The  microstructure  and  the  composition  of  the  catalyst  before  and  after  testing  were  analyzed  by 
scanning  electron  microscopy  (SEM,  Zeiss  FESEM)  combined  with  an  energy  dispersive  X-ray 
spectroscopy  (EDS).  The  catalytic  oxidation  of  carbon  monoxide  was  carried  out  in  a  horizontal 
flow  reactor,  illustrated  in  Fig.  1.  The  catalyst,  20g  MC/AI2O3  (approximately  0.2g  MC),  was 
loaded  to  the  center  of  the  alumina  tube  and  packed  tightly  with  two  porous  ceramic  pellets.  A 
mixuture  of  CO  and  air  was  fed  in  from  one  side  of  the  reactor  and  collected  from  the  other  side. 
The  flow  rate  of  each  gas  was  controlled  by  mass  flow  controllers  (Alicat  Scientific,  MC- 
50SCCM-D/5M),  while  soak  film  flow  meter  was  used  for  calibrations.  The  flow  rate  of  CO  was 
controlled  at  6.76  seem,  while  the  total  flow  of  approximately  100  seem.  The  actual  total  flow 
rate  of  inlet  gas  (feed  gas)  and  outlet  gas  (product  gas)  was  measured  by  the  soak  film  flow  meter 
and  used  for  calculating  the  CO  conversion.  The  outlet  gas  was  analyzed  by  a  CO2  detector 
(Vernier,  CO2  sensor)  to  identify  the  content  of  the  CO2.  The  concentration  of  CO2  was  tested 
every  10  °C  from  300-600  °C.  At  each  temperature  point,  the  system  was  held  at  least  for  30 
min  before  measuring  the  CO2  concentration.  The  CO2  content  is  then  used  to  determin  the  CO 
conversion  as  a  function  of  temperature. 

For  the  reaction  of  CO  oxidation,  CO  -1-  V2  O2  ^  CO2,  the  amount  of  CO  consumed  should  equal 
to  that  of  CO2  produced.  Therefore,  the  CO  conversion  as  percentage  can  be  calculated  by: 

Conversion  =F’co2/Fco=(  F’totaixCco2)/Fco , 

where  F’totai  is  outlet  gas  total  flow  rate,  F’co2  is  CO2  yield  flow  rate  in  outlet  gas,  Cco2  is  the 
concentration  of  CO2  in  the  outlet  gas,  and  Fco  is  the  CO  feed  flow  rate. 


Furnace 


Alumina  tube 


Figure  1:  Illustration  of  the  tubular  reactor  for  the  catalyst  testing 
D2.3:  DFT  Modeling 

All  structures  were  optimized  at  the  B3EYP  [29-30]/6-31G(d,p)  [31-34]  level  using  the  Gaussian 
09  program[35].  All  stationary  points  on  the  potential  energy  surface  (PES)  were  then  verified  by 
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calculated  vibrational  frequencies  at  the  same  theoretical  level.  Therefore,  the  reactant  and 
product  structures  are  truly  local  minima  while  each  transition  state  is  at  the  first  order  saddle 
point.  The  reaction  pathway,  i.e.  the  connection  between  reactant,  TS  and  product,  was  verified 
by  the  intrinsic  reaction  coordinate  (IRC)  calculations. 

D3:  Results  and  discussion 

D3.1:  Physical  characterization 

Fig.  2  (a)  and  (b)  shows  the  SEM  images  of  the  MC/AI2O3  catalyst  before  and  after  test, 
respectively.  The  flake  like  crystal  is  the  MC  phase,  which  covers  almost  all  the  surface  of  the 
AI2O3  support  comparing  to  pure  AI2O3,  as  shown  in  Fig.  2(c).  To  further  confirm  the 
composition  of  the  crystal  phase,  EDS  was  involved  to  test  the  composition  of  the  red  square 
area  in  Fig.  2(a),  and  the  result  was  shown  in  Fig.  2(d).  From  Fig.  2(d),  Na,  K,  C  and  O  were 
observed,  indicating  that  the  crystal  phase  is  the  molten  carbonate.  However,  Fi  is  too  light  to  be 
detected  by  EDS.  Element  A1  is  from  the  AI2O3  support,  and  Au  is  from  the  Au  coating,  which 
had  to  be  added  for  the  SEM  imagining.  A  schematic  of  the  MC/AI2O3  catalyst  below  melting 
point  is  shown  in  Fig.  2(e),  while  Fig.  2(f)  for  it  above  the  melting  point.  The  catalyst  was 
weighted  before  and  after  test  and  the  weight  was  almost  the  same,  indicating  no  any  loss  of  the 
catalyst  during  the  reaction  process.  In  addition,  its  microstructure  remains  same  before  and  after 
the  test.  So  the  MC/AI2O3  catalyst  has  a  good  stability. 


(f) 

Liquid  Phase 
Molten  Carbonate 


Above  melting  point 


Crystal  Phase 
Carbonate 


Figure  2:  SEM  images  of  (a)  MC/AI2O3  before  test;  (b)MC/Al203  after  test;  (c)Al203  pellet 
support;  (d)  EDS  result  of  MC  on  the  AI2O3;  (e)  schematic  of  solid  phase  of  MC  on  the 
AI2O3;  (f)  schematic  of  MC  on  the  AI2O3  above  melting  point 

D3.2:  Catalytic  reaction  evaluation 

Fig.  3  shows  the  conversion  of  the  earbon  monoxide  to  earbon  dioxide  with  and  without  the 
presence  of  the  MC/AI2O3  catalyst.  It  is  found  that  the  CO  conversion  with  MC/AI2O3  is  higher 
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in  the  whole  temperature  range  of  300-600  °C,  showing  a  good  catalytic  performance  of 
MC/AI2O3  for  the  oxidation  of  CO.  Two  inflection  points  are  shown  in  the  conversion  curve  of 
MC/AI2O3.  The  first  inflection  point  is  at  about  450  °C,  which  is  the  melting  point  of  the  MC. 
Below  450  °C,  the  CO  conversion  on  solid  phase  MC  is  very  low,  but  still  much  higher  than  that 
with  no  MC  catalyst.  However,  the  CO  conversion  on  liquid  phase  increases  dramatically  with 
the  temperature  increase  from  400  °C  to  450  “C.  The  CO  conversion  with  MC/AI2O3  reaches  96% 
at  500°C,  while  only  25%  at  the  same  temperature  without  catalyst.  The  results  incidate  that  both 
solid  phase  and  liquid  phase  MCs  present  the  catalytic  activity  for  the  CO  oxidation,  but  the 
liquid  phase  has  a  much  higher  catalytic  activity.  The  second  inflection  point  is  at  500°C,  the 
conversion  curve  becomes  flat  and  stabilized  at  96-98%  due  to  the  fact  that  most  of  CO  has  been 
converted  to  CO2. 


Figure  3:  CO  conversion  vs  Temperatrue 
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Figure  4:  Arrhenius  plot  of  calibrated  CO  conversion 
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Fig.  4  shows  the  arrhenius  plot  of  the  CO  conversion  by  MC/AI2O3  as  a  function  of  1000/T  (K"'). 
Before  molten  carbonate  melting,  the  activation  energy  (Ea)  is  approximately  27.0  kcal/mol. 
After  melting  point,  the  Ea  value  is  decreased  to  19.0  kcal/mol.  The  large  difference  in  Ea 
implies  that  the  reaction  may  proceed  through  different  pathways  under  low  and  high 
temperature.  Above  500  °C,  the  gas  phase  reaction  of  CO  with  oxygen  becomes  significant  and  it 
is  extremely  hard  to  evalue  the  contribution  from  the  catalyst.  So  we  did  not  examine  the  activity 
energy  beyond  500  °C. 

D3.3:  Reaction  mechanism 

Assuming  that  the  reaction  will  proceed  through  a  Eangmuir-Hinshelwood  (E-H)  mechanism,  it 
can  then  be  described  by  Eqn.  (1-3).  Using  DET  methods,  those  three  steps  have  been 
investigated  to  help  understand  the  reaction  kinetics  and  mechanism.  The  enthalpy  change  is 
calculated  to  be  -0.3  kcal/mol  and  -36.4  kcal/mol  at  the  B3EYP/6-31G(d)  level  for  Eqn.  (1)  and 
(2),  respectively.  This  indicates  that  both  oxygen  dissociation  and  CO  adsorption  are  exothermic 
and  favored  by  chemcial  equilibrium,  but  the  step  of  Eqn.  (2)  is  much  more  favored,  i.e.  K2  » 
Ki.  However,  the  step  of  Eqn.  (1)  has  an  energy  barrer  of  about  23.0  kcal/mol  calculated  by  DET 
in  pure  lithium  carbonate  [28]. 

+(C03)"-(MC)  ^^(C0,)"-(MC)  (1) 

CO(g)  +  (C03)"-(MC)  (2) 

(C^O^l'-fMC)  +  (C04)'-(MC)  - >2C03)"-(MC)  +CO3  (g)  (3) 

Eor  the  surface  reactioin  of  Eqn.  (3),  DET  calculations  were  performed  using  a  small  cluster  of 
Ei2(C204)-Ei2(C04)  as  shown  in  Eig.  5.  The  overall  reaction  path  and  potential  energy  surface 
(PES)  are  plotted  in  Eig.  6.  At  the  initial  state,  both  CO  and  O2  are  brought  to  the  surface  of 
carbonate  catalyst  as  the  molecular  form  of  (C204)^'  and  (C04)^'  and  the  oxygen  and  CO  are 
pointing  to  each  other,  highlighted  by  red  and  blue  circule,  respectively.  Once  the  oxygen  atom 
leaves  its  carbonate  carrier,  it  will  approach  CO  of  (C204)^'  and  this  leads  to  a  transition  state 

(TSl).  At  TSl,  the  0-0  bond  is  breaking,  while  the  C-0  bond  is  forming  as  highlighted  in  the 

orange  circle.  Next,  CO2  is  formed,  but  it  is  attached  to  a  carbonate  ion,  forming  an  intermediate 
state  of  pyrocarbonate  (C205)^'.  The  same  structure  was  identified  by  Huang  et  al.  and  reported 
in  Ref.  [36].  Through  another  transiton  state  (TS2),  CO2  is  released  and  finally  detached  from  the 
surface  of  MC.  Erom  Eig.  6,  the  energy  barrier  for  the  surface  reaction  of  Eqn.  (3)  is  calculated  to 
be  only  1 1.4  kcal/mol. 


(C;04)-(MO+(C04)^(MO 
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Figure  5:  Structureal  evolution  of  the  reaction  process  calculated  by  DFT 
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Figure  6:  Reaction  path  and  potential  energy  surface 

From  classical  Langmuir-Hinshelwood  mechanism,  the  reaction  rate  can  be  described  by: 


(4). 


where  r  is  the  reaction  rate,  k  for  the  rate  constant  of  Eqn.  (3),  Ki  and  K2  for  the  equilibrium 
constant  of  Eqn.  (1)  and  (2).  Erom  the  DET  results,  we  can  conclude  K^CQ^  then  the 

reaction  rate  will  be  approximated  as: 


il  +  K,C,o) 


(5). 


If  the  concentration  of  CO  is  extremely  low,  then  « 1  and  the  reaction  rate  is 

approximated  as 

r  =  kK,K,C'§CcoCl^.-  (6). 


The  reaction  is  the  first  order  with  respect  to  CO,  while  half  order  for  O2.  Increase  of  the 
concentration  of  reactant,  both  CO  and  O2,  will  increase  the  reaction  rate.  On  the  other  hand,  if 
the  CO  concentration  is  very  high,  then  the  reaction  rate  is  described  by  the 

equation  of 


kC 


K  C 


CO 


^  K  C 

^^2'^  CO 


(7). 


In  this  case,  the  CO  molecules  will  cover  the  catalyst  surface  and  will  inhibit  the  reaction.  This 
should  make  the  reactor  more  difficult  and  experience  higher  overral  energy  barrier.  In  the 
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current  experiment  setting,  the  effect  of  the  concentration  of  CO  was  not  investigated  yet  and 
will  be  reported  in  the  future. 

In  Fig.  6,  the  reaction  path  is  described  by  the  PES  as  a  function  of  reaction  coordinate.  The 
green  part  of  oxygen  dissociation  in  molten  carbonate  is  from  Ref.  [28].  The  reaction  is  initiated 
by  oxygen  dissociation,  and  then  undergoes  with  CO  adsorption,  surface  reaction  of  absorbed 
species,  and  finally  desorption  of  CO2.  The  rate  determing  step  is  oxygen  dissociative  adsorption 
with  an  energy  barrier  of  23.0  kcal/mol,  while  the  surface  reaction  is  relatively  easy  with 
calculated  activation  energy  of  11.4  kcal/mol.  The  actual  activation  energy  through  Arrheniuns 
plot  method  (Fig.  4)  is  measured  to  be  19.0  kcal/mol.  This  agrees  with  the  DFT  results.  It  should 
be  pointed  out  that  the  contribution  of  entropy  was  not  considered  in  the  DFT  calculations,  but  it 
should  have  some  contributions  to  ther  eaction  kinetics  in  such  a  high  temperature  reaction 
system. 

D4:  Conclusions 

In  the  current  study,  we  have  examined  molten  carbonate  as  a  non-metal  catalyst  for  the  CO 
oxidation  at  the  temperature  of  300-600  °C  for  the  first  time.  The  conversion  is  largely  inhanced 
by  melting  at  about  450  °C  and  increased  to  96%  at  500  °C.  The  reaction  mechanim  is  described 
by  a  Fangmuir-Hinshelwood  model  including  four  steps:  (1)  dissociation  adsorption  of  oxygen, 

(2)  adsorption  of  CO,  (3)  surface  reaction,  and  (4)  desorption  of  CO2.  DFT  modeling  reveals  the 
formation  of  (C204)^'  and  (C04)^'  as  the  intermediate  species.  It  also  shows  that  the  first  two 
steps  are  exothermic  and  preferred  by  chemical  equilibrium,  but  the  oxygen  dissociation  has 
higher  energy  barrer  than  the  surface  reaction  itself  and  become  the  rate  determing  step  in  the 
whole  process.  The  calculated  energy  barrier  for  oxygen  dissociation  to  form  C04^'  is  23.0 
kcal/mol  from  Ref.  [28],  which  is  in  a  good  agreement  with  the  measured  overral  activation 
energy  of  19.0  kcal/mol.  Further  analysis  of  the  reaction  kinetics  indicates  that  the  reaction  is 
affected  by  the  CO  concentration.  With  low  CO  concentration,  the  reaction  is  the  first  order  with 
respect  to  CO  and  half  order  to  O2.  From  this  study,  it  has  been  proven  that  molten  carbonate  can 
serve  as  an  efficient  catalyst  for  CO  oxidation  and  potentially  for  other  oxidation  reactions  at  the 
temperature  range  of  450-600  °C.  More  studies  are  demanded  to  further  deliberate  the  reaction 
mechanism  and  explore  more  potential  industrial  applications. 
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PART  E:  CO  OXIDATION  BY  ATOMIC  AND  MOLECULAR  OXYGEN  ON  AG(IOO) 

El:  Introduction 

Carbon  monoxide  (CO)  is  one  of  the  toxic  and  harmful  gases  in  the  air  and  also  a  poisoning  gas 
in  many  industrial  catalytic  processes.  Its  oxidation  by  selected  catalysts  has  been  of 
considerable  interest  because  of  its  significance  in  air  cleaning  and  gas  purification.  For  this 
purpose,  a  variety  of  catalysts  have  been  investigated.  Noble  metals  are  of  high  efficiency  at  low 
temperatures  but  their  applications  are  limited  by  high  costs  and  limited  amount  of  resources. 
Therefore,  more  research  has  been  turned  to  non-noble  metal-based  catalysts,  such  as  binary 
oxides  of  Co304/Ce02  CuO/CeOi  and  ternary  oxide  of  CuO/CosO/i/CeOi  (CCC)  Such 
oxides  above  have  shown  high  activity  for  the  CO  oxidation  in  simulated  exhaust  conditions. 
Additionally,  great  attention  has  recently  been  paid  to  Pt-based  intermetallic  compounds  For 
example,  the  Pt/Fe  bimetallic  nano-alloy  is  favorable  for  CO  oxidation,  over  which  CO  can  be 
completely  removed  around  the  room  temperature.  Besides  metal  oxide  and  bimetallic  catalysts, 
new  efforts  are  also  being  made  to  search  for  better  supporting  materials  such  as  monolayer 
M0S2  graphyne  graphene  graphene  oxide  and  single  layer  hexagonal  BN  (h-BN) 
They  can  not  only  effectively  anchor  the  metal  atoms,  but  also  can  boost  the  catalytic 
performance  due  to  their  large  specific  surface  areas  and  novel  electronic  and  thermal  properties. 
This  opens  a  new  direction  in  designing  future  catalysts  for  the  CO  oxidation. 

Compared  to  catalysts  mentioned  above,  silver  (Ag)  as  a  relatively  inexpensive  metal  is  an 
important  candidate  for  the  CO  oxidation  reaction  in  present  industries.  Ag-based  catalysts  have 
been  extensively  investigated  in  recent  years  due  to  their  catalytic  properties  for  CO  oxidation 
The  results  in  Ref.  [24]  show  that  porous  Ag-Fe203  nanocomposite  has  a  good  reproducibility 
and  stability  in  CO  oxidation.  Zhang  et  al.  reported  that  an  Ag/SBA-15  catalyst  with  low  loading 
of  1.42  wt%  can  achieve  98%  conversion  of  CO  at  20  °C  More  recently,  it  was  shown  that 
well-dispersed  spherical  silver  nanoparticles  deposited  on  carbon  nanotube  can  convert  CO  to 
CO2  approximately  100%  at  220  °C  On  the  theoretical  side,  the  CO  oxidation  was  also 
examined  on  silver  or  silver-based  alloys  using  cluster  models  in  combination  with  density 
functional  theory  (DFT)  methods  It  is  surprising  that  simple  silver  surfaces  such  as  100,  110 
and  111  have  not  been  considered  yet,  which  are  of  considerable  importance  for  a  better 
understanding  of  the  reaction  kinetics  and  mechanism.  In  this  report,  a  systematic  DFT  study  of 
CO  oxidation  by  molecular  and  atomic  oxygen  on  the  surface  of  Ag(lOO)  is  presented.  The  goal 
is  to  obtain  the  mechanism  for  CO  oxidation  at  the  atomic  level.  Hence,  the  reaction  pathways, 
energetics,  geometrical  and  electronic  structures  (density  of  state  and  molecular  orbitals)  were 
particularly  computed  and  discussed. 

E2:  Computational  details  and  models 

All  calculations  were  performed  using  the  suite  program  of  Gaussian09  The  pw91pw91 
hybrid  exchange-correlation  functional  was  used,  which  has  shown  good  reliability  and 
accuracy  in  previous  reports  For  silver,  the  effective  core  potential  (ECP)  of  Lanl2 
combined  with  the  DZ  basis  set  was  adopted  While  for  O  and  C  atoms,  the  full  electron 
basis  set  of  6-31G(d)  was  employed  All  stationary  points  on  the  potential  energy  surface 
(PES)  were  verified  by  calculated  vibrational  frequencies  at  the  same  theoretical  level  as  used  in 
the  geometry  optimization.  To  locate  a  transition  state  (TS),  a  relaxed  PES  scan  was  carried  out 
first  to  find  the  structure  with  the  highest  energy  on  the  reaction  path.  It  was  then  optimized  with 
full  relaxation  using  the  Bemy  algorithm  and  confirmed  to  be  a  saddle  point  with  only  one 


42 


PROPOSAL  NUMBER:  62940-CH-REP,  FINAL  REPORT  2013-2017 


imagery  frequency.  Finally,  the  reaction  path  was  determined  by  the  intrinsic  reaction  coordinate 
(IRC)  calculations,  which  verifies  the  connection  from  the  reactant,  TS  to  product  on  PES. 

In  the  current  work,  only  the  low  index  surface  of  Ag(lOO)  was  considered  as  the  model  for 

silver  surface,  which  was  represented  by  a  three-layer  silver  cluster  with  a  fixed  Ag-Ag  distance 

0 

of  2.89  A,  the  value  from  experimental  bulk  lattice.  As  displayed  in  figure  1,  the  cluster  models 
of  Ag29  and  Agss  represent  the  Ag(lOO)  surface  and  are  used  for  CO  oxidation  by  atomic  oxygen 
and  molecular  oxygen,  respectively.  Three  adsorption  sites  were  considered:  on  the  top  of  silver 
atom  (atop  site),  on  the  top  of  Ag-Ag  bond  (bridge  site),  and  on  the  top  of  the  center  of 
tetragonal  silver  ring  (hollow  site).  The  adsorption  energy  is  defined  by  the  formula  of 

^ ads  ^ adsorbate~^  ^ ^ adsorbatd A^ilOO) 

where  represents  the  total  energy  of  O,  triplet  O2,  and  CO,  respectively.  EAgcioo) 

^adsorbatiAg(wo)  ^^c  total  cncrgy  of  Ag(lOO)  without  and  with  the  adsorbate.  By  this  definition, 
a  positive  value  of  adsorption  energy  represents  an  exothermic  adsorption  process. 

To  test  the  reliability  of  pw91pw91,  a  comparative  study  of  the  adsorption  of  CO  on  Ag(lOO) 
were  performed  using  b3pw91  and  bSlyp  The  results  are  presented  in  subsection  3.1  with 
adsorption  energies  summarized  in  Table  1.  From  table  1,  it  is  found  that  the  values  by 
pw91pw91  are  in  the  best  agreement  with  those  from  experiments  and  other  theoretical 
calculations.  In  addition,  we  did  some  testing  calculations  on  the  CO  adsorption  using  long -rang 
corrected  B3LYP  functionals  (CAM-B3LYP,  GD2/B3LYP  and  GD3/B3LYP).  The  results  are 
listed  in  Table  2.  It  was  found  that  GD2  and  GD3  can  largely  improve  the  adsorption  energy,  but 
have  only  small  effects  on  the  natural  charge  and  adsorption  structures.  So  those  testing  results 
have  confirmed  that  pw91pw91  is  the  best  choice  for  the  current  molecular  system. 


1-fold 


Fig.  1  Cluster  models  (Ag29  and  Agss)  for  Ag(lOO)  and  three  typical  adsorption  sites  of  atop, 
bridge,  and  hollow.  The  model  is  designed  by  a  three-layer  silver  cluster  with  fixed  Ag-Ag 

o 

distance  of  2.89  A.  Ag29  and  Agss  represent  the  surface  of  Ag(lOO)  for  the  CO  oxidation  by 
atomic  and  molecular  oxygen,  respectively. 


Tablel.  The  calculated  C-O  bond  length  (dc-o),  CO-surface  distance  of  (dc-Ag),  C-O 
stretching  frequency  (coc-o),  CO-Ag  vibrational  frequency  (cuco-Ag),  and  adsorption  energy 
(Eads)  for  the  CO  adsorption  on  Ag(lOO)  together  with  previous  experimental  and 
theoretical  values. _ 

dc-o  (A)  dc-Ag  (A)  coc-o  (cm'^)  cuco-Ag  (cm'^)  Eads 

(kcal/mol) 
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pw91pw91  1.152  2.232  2055.5  156.5  8.4 


b3pw91  1.138 

2.312  2176.0 

123.6 

3.7 

b31yp  1.137 

2.469  2181.5 

83.4 

2.7 

theoretical  1.150 

52  2.230^2  2050^2 

203 

52 

experimental  - 

2137.052 

- 

6.554 

Table2.  The  calculated  natural  charge  on  C  atom  (Qc),  on  O  atom  (Qo),  sum  of  CO 
molecule  (Qco),  C-O  bond  length  (dc-o),  CO-surface  distance  of  (dc-Ag),  and  adsorption 
energies  of  CO  (Eads)  on  silver  (Agss)  obtained  from  PW91PW91,  B3LYP,  CAM-B3LYP, 
GD2/B3LYP,  and  GD3/B3LYP,  respectively. _ 


Qc  (e) 

Qo  (e) 

Qco  (e) 

dc-o(A) 

dc-Ag(A) 

Eads  (kcal/mol) 

PW91PW91 

0.559 

-0.374 

0.185 

1.152 

2.232 

8.4 

B3LYP  0. 

,638  -1 

0.423  0 

.215  1, 

.137  2.‘ 

469 

2. 

1 

CAM-B3LYP 

0.648 

-0.421 

0.227 

1.132 

2.409 

4.6 

GD2/B3LYP 

0.613 

-0.409 

0.204 

1.138 

2.301 

8.0 

GD3/B3LYP 

0.620 

-0.412 

0.208 

1.137 

2.334 

8.5 

E3:  Results  and  Discussions 

E3.1:  Adsorptions  of  single  O  atom,  O2  and  CO  molecules  on  Ag(lOO) 

We  first  examined  the  adsorption  of  single  O  atom,  O2  and  CO  molecules  on  the  clean  surface  of 
Ag(lOO),  which  is  the  prerequisite  for  the  CO  oxidation.  Various  adsorption  configurations  have 
been  considered  in  order  to  find  the  most  stable  one  for  each  adsorbate.  Fig.  2  shows  the 
optimized  structures  with  the  highest  adsorption  energy  for  single  O,  O2  and  CO,  respectively. 
For  atomic  oxygen,  it  prefers  the  4-fold  hollow  site  with  an  adsorption  energy  of  76.6  kcal/mol. 
While  for  molecular  oxygen,  the  O2  molecule  is  parallel  to  the  silver  substrate  and  adsorbed  on 
the  hollow  site  with  only  small  adsorption  energy  of  2.7  kcal/mol.  The  distance  between  oxygen 
and  silver  surface  (do-Surf)  is  1.037  A  and  1.724  A  for  atomic  and  molecular  oxygen,  respectively. 
From  the  calculated  adsorption  energy  and  do-surf,  we  conclude  that  the  adsorption  of  O  atom  on 
Ag(lOO)  is  stronger  than  that  of  O2  molecule.  In  addition,  the  natural  population  analysis  (NPA) 
indicates  that  the  sum  of  natural  charges  is  -0.943  e  and  -0.784  e  for  atomic  and  molecular 
oxygen,  respectively.  The  negative  charge  implies  a  significant  charge  transfer  from  silver 
substrate  to  the  O  atom  or  O2  molecule,  showing  a  strong  chemical  bonding  between  silver  and 
oxygen.  So  the  adsorption  of  either  atomic  or  molecular  oxygen  is  of  chemical  type  on  Ag(lOO). 

o 

Furthermore,  after  O2  adsorption,  the  0-0  bond  length  is  elongated  to  1.397  A,  which  is 
significantly  larger  than  that  of  the  free  O2  molecule,  1.228  A.  Compared  with  the  vibrational 
frequency  of  free  O2  molecule  (1574.6  cm'*),  the  0-0  stretching  frequency  of  838.7  cm'*  is 
drastically  lowered  by  the  bonding  to  the  silver  surface.  In  the  view  of  the  elongated  0-0  bond 
and  very  low  vibrational  frequency,  it  is  clear  that  the  0-0  bond  is  weakened  and  the  O2 
molecule  is  likely  turned  into  superoxide  on  silver.  This  result  is  very  similar  to  the  dioxygen  on 
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Ag  (110)  reported  by  Backx  et  al.  The  weakening  of  the  0-0  bond  makes  dioxygen  more 
reactive,  which  will  facilitate  the  oxidation  of  CO  molecule  on  the  silver  surface. 

Different  from  oxygen,  CO  is  found  to  be  the  most  stable  at  the  1-fold  atop  site  in  an  upright 
position  with  the  carbon  atom  pointing  towards  the  silver  surface  (C-down  adsorption  model), 
which  agrees  with  previous  studies  using  the  configuration  interaction  (Cl)  method  The 
corresponding  adsorption  energy  and  C-Ag  distance  is  8.4  kcal/mol  and  2.232  A,  respectively. 
The  present  adsorption  energy  is  very  close  to  the  experimental  value  of  6.5  kcal/mol  obtained 
using  thermal  desorption  spectroscopy  (TDS)  on  Ag  (111)  surface  but  it  is  slightly  higher  than 
4.4  kcal/mol  at  the  Cl  level  on  Ag(lOO)  and  ~3.2-3.7  kcal/mol  at  the  DFT-GGA-PW91  level 
on  Ag(l  1 1).  The  net  natural  charge  transfer  from  CO  to  Ag(lOO)  is  found  to  be  0. 185  e  (C:  0.559 
e;  O:  -0.374  e),  which  is  slightly  larger  than  0.123  e  obtained  from  Cl  calculations  and  mainly 
from  the  C  lone  pair  orbital  to  silver  The  calculated  molecular  orbital  in  Fig.  3  confirms  that 
CO  adsorption  on  silver  is  of  chemical  type.  Upon  adsorption,  the  C-0  bond  length  (dc-o)  is 
elongated  due  to  the  interaction  between  C  2p  and  Ag  Ad  orbitals.  The  dc-o  is  calculated  to  be 

o  o 

1.152  A,  which  is  slightly  larger  than  that  of  free  CO  molecule,  1.148  A,  but  in  line  with  the 
1.150  A  calculated  using  the  projector- augmented-wave  (PAW)  method  .  Moreover,  the 
vibrational  frequency  of  C-0  is  calculated  to  be  2055.5  cm'\  which  is  smaller  than  that  of  free 
CO  molecule  (2128.1  cm'*)  and  experimental  value  of  2137.0  cm'* 


Fig.  2  The  optimized  adsorption  structure  with  the  lowest  energy  for  single  O  atom  (a),  O2 
molecule  (h),  and  CO  molecule  (c)  on  Ag(lOO). 


Fig.  3  Selected  molecular  orbital  of  CO  adsorption  on  Ag(lOO). 
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E3.2:  CO  oxidation  by  molecular  oxygen  on  Ag(lOO) 

As  we  know,  there  are  two  well-established  mechanisms  for  the  CO  oxidation  by  O2  molecule: 
the  Langmuir- Hinshelwood  (L-H)  13,57-58  Eley-Rideal  (E-R)  mechanism  i*’39-6i 

E-H  mechanism,  the  CO  and  O2  molecules  are  co-adsorbed  on  the  surface  of  silver  before 
reaction,  a  peroxo-type  OOCO  complex  is  proposed  to  be  an  intermediate  state  (MS).  However, 
in  the  E-R  mechanism,  the  CO  molecule  directly  reacts  with  the  pre-adsorbed  O2  molecule, 
producing  a  carbonate-like  CO3  complex  as  the  MS.  In  the  present  study,  the  calculated 
adsorption  energies  for  both  CO  (8.4  kcal/mol)  and  O2  (2.7  kcal/mol)  on  Ag(lOO)  are  positive. 
So  the  E-H  mechanism  is  expected  to  prevail  over  the  E-R  mechanism  for  CO  oxidation  by  O2 
molecule  on  the  silver  surface.  Erom  previous  studies  the  complete  CO  oxidation  includes  two 
steps:  the  E-H  mechanism  as  a  starting  step  with  the  reaction  of  C0-i-02^00C0^C02-i-0,  then 
following  the  E-R  mechanism  with  the  reaction  of  C0+0^C02.  To  achieve  a  comprehensive 
and  better  understanding,  both  mechanisms  of  CO  oxidation  are  discussed  below. 

Eirst,  we  considered  the  E-H  mechanism  of  CO  oxidation.  Eigure  4  depicts  the  structures  of 
initial  state  (IS),  transition  state  (TS),  and  final  state  (ES)  together  with  the  energy  barrier  of  CO 
oxidation  on  Ag(lOO)  with  O2  molecule  present.  As  seen  from  figure  4,  after  CO  and  O2  co- 

o 

adsorbed  on  silver,  CO  is  nearly  parallel  to  O2  with  the  C-02  distance  of  1.356  A  and  the  02-C- 
01  angle  of  119.3°.  At  this  phase,  the  bond  length  of  02-03  (1.461  A)  is  larger  than  that  of  free 
oxygen  molecule,  indicating  that  the  adsorbed  oxygen  molecule  is  activated  and  ready  to 
participate  in  the  CO  oxidation  reaction.  The  01-C  bond  length  is  1.229  A,  slightly  larger  than 
that  in  the  free  CO  molecule. 

As  displayed  in  figure  4,  the  CO  oxidation  process  is  initiated  by  the  stretching  of  oxygen 
molecule  (02-03)  and  leads  to  a  TS  on  PES.  At  TS  (only  imaginary  frequency  is  /465.6  cm'^), 
the  02-03  distance  is  elongated  to  1.793  A  and  the  C-02  distance  is  shortened  to  1.280  A, 
showing  a  simultaneous  bond  forming  and  breaking  process.  At  the  end  of  reaction,  C-01  and  C- 
02  bond  lengths  are  shortened  to  the  equal  value  of  1.179  A,  which  is  the  same  as  the  bond 
length  in  free  CO2.  The  01 -C-02  angle  becomes  180.0°  and  a  linear  CO2  molecule  is  formed. 
The  calculated  energy  barrier  is  only  6.9  kcal/mol,  indicating  a  favorable  reaction  of  CO 
oxidation  on  silver.  In  addition,  it  is  noticeable  that  the  bond  of  01-C  is  first  elongated  from 

00  o 

initial  state  of  1.229  A  to  transition  state  of  1.244  A,  and  then  shortened  to  1.179  A  at  final  state. 

o  o 

The  distance  between  03  and  Ag  surface  decreases  from  1.8  A  of  initial  state  to  1.6  A  of 
transition  state  to  1.1  A  of  final  state.  The  variation  of  such  bond  distances  sufficiently  indicates 
that  the  CO  oxidation  is  a  synergistic  process  with  CO,  O2  and  silver  surface  all  actively 
involved. 

E3.3:  CO  oxidation  by  atomic  oxygen  on  Ag(lOO) 

In  this  section,  we  will  focus  on  the  E-R  mechanism  of  CO  oxidation.  As  discussed  above,  the 
final  state  of  CO  oxidation  by  O2  molecule  leads  to  a  detached  CO2  and  another  chemisorbed  O 
atom  on  the  silver  surface.  Next,  the  CO  oxidation  by  atomic  oxygen  on  the  silver  surface  will  be 
examined.  Eigure  5  shows  the  local  structures  of  beginning,  transitioning,  and  ending  states  for 
the  CO  oxidation  by  atomic  oxygen  on  the  Ag(lOO)  surface.  Eirst,  CO  is  adsorbed  on  Ag(lOO) 
with  an  O  atom  pre-adsorbed.  The  calculated  adsorption  energy  of  CO  is  6.2  kcal/mol  by 
Eads  =  Eco(gas)  +  ^’o/Agdoo)  “  ^’co-o/ Ag(ioo)  >  whlch  Is  sllghtly  lower  than  that  calculated  on  the 

clean  Ag(lOO).  However,  both  values  are  in  fairly  good  agreement  with  the  value  of  6.5  kcal/mol 
observed  by  experiment  for  CO  molecule  adsorbed  on  the  clean  surface  of  Ag(lOO),  indicating 
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small  effect  of  pre-adsorbed  oxygen  on  the  adsorption  of  CO.  In  figure  5,  at  the  beginning  state, 
the  CO  molecule  is  weakly  adsorbed  on  the  top  of  Ag2  with  a  tilt  angle  of  25.6°  and  C-02 
distance  of  2.796  A.  The  C-01  bond  length  of  1.149  A  is  very  close  to  the  1.148  A  for  free  CO 
molecule  and  slightly  shorter  than  1.152  A  for  CO  adsorption  on  the  clean  silver  surface. 

o 

Compared  with  the  adsorption  on  the  clean  silver  surface,  the  C-Ag2  distance  of  2.570  A  is 
larger  than  that  of  2.243  A,  while  the  O2-Ag(100)  distance  of  1.003  A  is  slightly  short  than  that 
of  1.085  A. 


Fig.  4  Optimized  structures  and  activation  energy  for  CO  oxidation  by  molecular  oxygen 
on  Ag(lOO). 

As  the  CO  molecule  (COl)  descends  and  approaches  the  surface,  the  distance  of  C-Ag2 
decreases.  At  the  same  time,  the  CO  molecule  rotates  itself  and  the  02  atom  accordingly  raises 
itself.  As  a  result,  the  01 -C-02  angle  increases  and  the  C-02  distance  decreases  quickly.  This 
leads  to  a  TS  located  on  the  PES  with  a  single  imaginary  frequency  of  /H8  cm"'.  At  TS,  the  C- 
02  distance  is  significantly  shortened  to  2.082  A  from  2.796  A,  while  the  C-01  bond  length  is 
enlarged  to  1.162  A,  and  the  Ol-C-02  angle  is  increased  to  113.7°  from  the  103.8°  at  the 
beginning  state.  Simultaneously,  the  C-Ag2  distance  decreased  to  2.382  A  while  the  02-surface 

o 

distance  increased  to  1.158  A.  During  this  process,  the  system  experiences  an  energy  barrier  of 
only  2.1  kcal/mol.  Such  a  small  conversion  barrier  implies  that  it  is  very  easy  for  CO  oxidation 
by  atomic  oxygen  on  the  silver  surface  disregarding  the  effect  of  entropy.  Finally,  the  01 -C-02 
angle  is  increased  to  179.9°,  and  the  bond  lengths  of  01-C  and  C-02  equal  to  1.179  A,  nearly 

o 

identical  to  those  in  the  free  CO2  molecule  (bond  angle:  180.0°;  bond  length:  1.179  A).  This 
indicates  the  formation  of  CO2  and  a  complete  oxidation  of  CO.  The  PES  calculated  by  IRC 
shows  the  connection  from  the  starting  structure,  TS,  to  the  final  structure,  providing  a  full 
picture  of  CO  oxidation  to  CO2  on  the  Ag(lOO)  surface  with  O  atom  pre-adsorbed. 
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Fig.  5  Optimized  structures  and  activation  energy  for  CO  oxidation  by  atomic  oxygen  on 
Ag(lOO). 


E3.4:  CO  oxidation  mechanism 

E3.4.1  CO  oxidation  mechanism  by  molecular  oxygen 

To  further  understand  the  mechanism  of  CO  oxidation  by  molecular  and  atomic  oxygen  on  the 
silver  surface,  we  plotted  the  broadened  density  of  state  (DOS)  and  overlap  DOSs  as  shown  in 
figure  6  and  7,  respectively.  Here,  the  electronic  DOS  was  obtained  by  broadening  the  electronic 
energy  levels  with  a  Gaussian  function.  The  full  width  at  half  maximum  (FWHM)  and  the  scale 
ratio  for  DOS  were  set  to  be  0.05  and  0.1,  respectively.  In  addition,  the  isosurfaces  of  selected 
molecular  orbitals  (MOs)  in  the  process  of  CO  oxidation  are  also  displayed  in  figure  6  and  7  with 
an  isovalue  of  0.02. 

For  CO  oxidation  by  oxygen  molecule  on  silver,  figure  6  shows  the  broadened  DOS  and  overlap 
DOSs  of  each  atom.  01  (red  line),  02  (blue  line)  and  03  (dark  teal)  represented  the  DOS  of  the 
O  atom  in  the  adsorbed  CO  molecule,  one  O  atom  in  the  adsorbed  O2  molecule  finally  bonded  to 
the  C  atom,  and  another  O  atom  in  the  adsorbed  O2  molecule  finally  chemisorbed  on  the  silver 
surface,  respectively.  OPl  (dark  yellow  line),  OP2  (dark  blue  line)  and  OPS  (brown  line) 
represented  the  overlap  DOSs  of  C-01,  C-02  and  02-03,  respectively.  The  black  line  and  pink 
line  were  DOSs  of  C  atom  and  Ag  atom,  respectively.  Be  aware  that  the  overlap  of  DOS  (OPn)  is 
greater  or  less  than  zero  denotes  the  corresponding  MO  is  favorable  or  unfavorable  for  forming 
chemical  bond,  respectively.  Different  from  the  CO  oxidation  by  atomic  oxygen,  the  process  of 
CO  oxidation  by  molecular  oxygen  involves  the  dissociation  of  oxygen  molecule  on  the  silver 
surface.  In  other  words,  the  0-0  bond  needs  to  be  broken  as  the  CO  molecule  being  oxidized.  As 
seen  from  DOSs  and  overlap  DOSs  in  figure  6,  at  the  reactant  phase,  there  are  overlaps  of  DOS 
between  C  and  01,  C  and  02,  02  and  03  at  the  energy  range  of  -13  eV  -  -10  eV,  which  is  also 
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evidenced  from  OPl,  OP2,  and  OPS.  So  the  orbitals  are  conductive  for  bonding  in  this  energy 
range,  which  is  also  observed  in  the  MOs.  For  the  instance  of  MO  on  the  left,  there  is  a  7i-like 
bond  between  C,  01  and  02.  Furthermore,  there  are  significant  hybrid  orbitals  between  s  and  p 
orbitals  of  02  and  03  atom,  resulting  in  a  a  bond  between  02  and  03.  Likewise,  there  are  5,  p 
and  d  hybrid  orbitals  between  C,  03  and  Ag,  leading  to  partial  a-bonding  between  03  and  Ag,  C 
and  Ag,  respectively.  For  the  MO  on  the  right,  it  is  mainly  of  a  big  7i-bond  between  01,  C,  02, 
and  03,  mixed  with  partial  a-bond  of  03-Ag,  and  C-Ag,  respectively. 

At  the  TS,  the  overlap  of  DOS  between  C  and  02  is  increased  while  that  between  02  and  03 
decreased,  implying  that  the  bonding  between  C  and  02  is  enhanced  and  that  between  02  and 
03  is  weakened.  On  the  other  hand,  the  interaction  between  C  and  Ag  is  largely  decreased  as  the 
C  atom  is  leaving  the  silver  surface  and  the  Ol-C-02  bond  angle  being  increased.  On  the 
contrary,  the  interaction  between  03  and  Ag  is  largely  enhanced  with  the  03  atom  getting  closer 
to  the  surface  and  the  bond  length  of  02-03  being  elongated.  The  results  from  DOS  analysis  are 
also  supported  by  the  selected  MOs  at  TS. 

Finally,  the  0-0  bond  of  oxygen  molecule  is  completely  broken  and  CO  oxidation  is  finished. 
From  the  DOSs  and  overlap  DOSs  of  product,  it  is  found  that  the  DOSs  of  01  and  02  are 
degenerate,  so  are  the  overlap  DOSs  of  OPl  and  0P2,  but  the  overlap  DOSs  of  0P3  is  almost 
zero  in  the  whole  considered  energy  range.  Correspondingly,  the  MO  (left)  at  -12.4  eV  forms  a  n 
bond  by  2p  orbitals  of  01,  C,  and  02  atom,  which  is  parallel  with  the  nodal  plane.  The  MO 
(right)  at  -12.3  eV  also  forms  a  n  bond,  which  is  perpendicular  to  the  nodal  plane. 


Fig.  6  The  density  of  state  (DOS),  overlap  DOS,  and  the  selected  molecular  orhitals  (MOs) 
related  to  the  bonding  of  01-C-02-03-Ag  at  reactant,  transition  state  (TS),  and  product 
for  CO  oxidation  hy  molecular  oxygen  on  Ag  (100),  respectively.  Ol,  02  and  03  represent 
the  O  atom  in  the  adsorbed  CO  molecule,  one  O  atom  in  the  adsorbed  O2  molecule  finally 
bonding  to  the  C  atom,  and  another  O  atom  in  the  adsorbed  O2  molecule  finally 
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chemisorbed  on  the  Ag  surface,  respectively.  OPl,  OP2  and  OP3  represent  the  overlap 
DOSs  of  C-Ol,  C-02  and  02-03,  respectively.  Red  and  green  parts  correspond  to  positive 
and  negative  regions  respectively,  and  the  iso  value  of  the  isosurface  is  set  to  be  0.02.  For 
DOS  of  product,  the  DOSs  of  01  and  02  atoms  are  degenerate,  so  are  the  overlap  DOSs  of 
OPl  and  OP2  are  degenerate,  but  the  overlap  DOSs  of  OP3  is  almost  zero  in  the  whole 
considered  energy  range. 


E3.4.2:  CO  oxidation  mechanism  by  atomic  oxygen 

For  CO  oxidation  by  atomic  oxygen,  as  shown  in  figure  7,  01  (red  line)  and  02  (blue  line)  are 
defined  as  the  DOSs  of  the  oxygen  in  the  adsorbed  CO  molecule  and  one  directly  adsorbed  on 
the  silver  surface,  respectively.  OPl  (pink  line)  and  0P2  (dark  yellow  line)  are  defined  as  the 
overlap  DOSs  of  C-01  and  C-02,  respectively.  At  the  beginning  of  CO  oxidation  (reactant),  the 
red  line  is  high,  showing  overlaps  with  the  black  line  (DOS  of  C)  in  the  region  of  -15  eV  to  -10 
eV.  However,  the  dark  teal  curve  (DOS  of  Ag)  and  the  blue  one  mainly  appear  in  the  region  of  - 
10  eV  to  HOMO.  The  pink  curve  is  greater  than  zero  in  the  region  of  -12.5  eV  to  -10  eV  under 
the  HOMO  energy  lever,  while  the  dark  yellow  curve  is  almost  zero  in  the  whole  considered 
energy  range.  Therefore,  we  can  conclude  that  2p  orbitals  of  01  and  C  have  made  significant 
contributions  to  the  bonding  between  C  and  01.  This  is  also  evidenced  in  the  corresponding 
MOs.  In  addition,  Ag  4d  and  02  2p  orbitals  have  significant  overlaps  in  the  band  of  -10  eV  to 
HOMO.  On  the  other  hand,  there  is  only  very  little  interactions  between  C  2p  and  02  2p  orbitals 
since  the  overlap  DOS  between  C  and  02  (OP2)  is  almost  zero  in  our  considered  energy  range. 
These  conclusions  are  then  confirmed  by  the  isosurfaces  of  corresponding  MOs.  For  example, 
the  MO  on  the  left  is  mainly  attributed  to  a  C-01  7i-bond,  while  that  on  the  right  is  a  7i-bond 
perpendicular  to  the  nodal  plane.  All  these  facts  reveal  that  the  C-01  remains  the  bonding 
characteristics  of  a  free  CO  molecule  and  is  positioned  for  bonding  with  02  atom  on  the  silver 
surface.  All  these  facts  support  an  E-R  reaction  mechanism  for  this  reaction  step.  It  is  necessary 
to  make  it  clear  that  the  effects  of  temperature  are  not  considered  here,  which  could  make  the 
free  energy  change  of  the  CO  adsorption  positive  due  to  the  negative  change  of  entropy  at  high 
temperatures.  The  actual  mechanism  of  the  CO  oxidation  by  atomic  oxygen  is  largely  affected  by 
the  temperature.  Therefore,  the  conclusion  of  an  E-R  mechanism  here  should  not  be  exaggerated 
and  further  examinations  by  molecular  dynamics  (MD)  based  methods  are  needed  to  clarify  the 
exact  reaction  mode. 

At  TS,  the  overlap  between  red  and  black  lines  is  reduced,  while  that  between  blue  and  black 
lines  is  increased  in  the  energy  window  of  [-11.5,  -8.5  eV].  This  is  supported  by  the  decreasing 
of  peak  intensity  of  OPl.  At  the  same  time,  the  OP2  intensity  appears  greater  than  zero  in  the 
same  energy  range.  These  evidences  indicate  that  the  bond  of  C-01  is  weakened  and  the  bond  of 
C-02  is  strengthened.  Meanwhile,  the  interaction  between  C  atom  and  Ag  surface  is  increased. 
Prom  the  MO  on  the  left,  as  the  CO  molecule  approaches  the  02  atom  and  Ag  surface,  a  a-bond 
is  formed  by  C  2p  and  02  2p.  The  overlap  between  Ag  4d  and  C  2p  is  also  significant.  The 
distinct  interactions  of  C-02  and  C-Ag  are  also  observed  in  the  MO  on  the  right.  These  bonding 
and  anti-bonding  characteristics  of  MOs  are  coincident  with  the  above  analysis  of  bond  length 
and  bond  angle  of  Ol-C-02.  At  the  ending  of  CO  oxidation  (product),  as  the  01,  C  and  02 
adjust  their  bond  length,  bond  angle  and  orientation,  the  CO2  molecule  is  finally  formed  and 
become  parallel  to  the  surface.  In  figure  6,  the  DOSs  of  01  and  02  atoms  are  degenerate,  so  are 
the  overlap  DOSs  of  OPl  and  0P2. 
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Fig.  7  The  density  of  state  (DOS),  overlap  DOS,  and  selected  molecular  orbitals  (MOs) 
related  to  the  bonding  of  01-C-02  at  reactant,  transition  state  (TS),  and  product  for  CO 
oxidation  by  atomic  oxygen  on  Ag  (100),  respectively.  01  and  02  represent  the  O  atom  in 
the  adsorbed  CO  molecule  and  the  O  atom  pre-adsorbed  on  silver,  respectively.  OPl  and 
OP2  represent  the  overlap  DOS  between  C  and  Ol,  and  between  C  and  02,  respectively. 
Red  and  green  parts  correspond  to  positive  and  negative  regions  respectively.  The  isovalue 
of  the  isosurface  was  set  to  be  0.02.  For  DOS  of  product,  the  DOSs  of  Ol  and  02  are 
degenerate,  and  so  are  the  overlap  DOSs  of  OPl  and  OP2. 


E3.4.3:  Formation  and  Dissociation  of  CO3  on  Ag(lOO) 

For  the  reaction  of  CO  oxidation  on  silver,  the  formation  of  surface  carbonate  species 
(C02(g)+0(ad)^C03(ad))  has  been  identified  Two  adsorption  structures  are  bi-coordinated 
perpendicular  and  tri-coordinated  parallel  with  calculated  adsorption  energy  of  91.2  and  89.1 
kcal/mol,  respectively.  For  the  dissociated  reaction  C03(ad)^0(ad)+C02(g),  the  energy  barrier  is 
calculated  to  be  19.8  and  26.3  kcal/mol  for  the  bi-coord  and  tri-coord  pathways,  respectively. 
However,  the  reaction  of  C03(ad)-i-C0(ad)^2C02(g)  has  an  energy  barrier  of  only  15.0  kcal/mol  in 
Ref.  [31].  In  the  same  study,  the  activation  energy  for  0(ad)-i-C0(ad)^C02(g)  was  calculated  to  be 
1.2  kcal/mol,  which  is  very  close  to  2.1  kcal/mol  from  the  current  study.  Based  on  those  facts,  we 
can  conclude  that  the  carbonates  are  fairly  stable  on  silver,  but  can  be  removed  by  further  CO 
dosing  or  dissociation.  The  reaction  pathway  of  CO  oxidation  through  carbonate  is  less  favored 
thermodynamically. 


E4:  Conclusions 
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In  summary,  we  have  studied  the  CO  oxidation  by  molecular  and  atomic  oxygen  on  the  silver 
surface  at  the  DFT-PW91PW91/[LANL2DZ,  6-31G(d)]  level.  The  Ag29  and  Ag38  cluster  models 
were  selected  to  represent  the  Ag(lOO)  surface.  The  calculated  results  show  that  all  adsorptions 
of  O,  O2  and  CO  molecules  on  Ag  (100)  are  of  chemical  type.  The  O2  molecule  is  partially 
activated  when  bonded  to  the  Ag  surface  and  ready  to  participate  in  the  reaction  of  CO  oxidation. 
The  CO  oxidation  by  O2  molecule  prefers  the  L-H  mechanism  and  it  involves  the  dissociation  of 
the  0-0  bond  and  the  formation  of  C-0  bond.  The  interactions  between  2s  and  2p  of  C  and  O, 
and  4d  of  Ag  orbitals  facilitate  the  CO  oxidation.  However,  the  CO  oxidation  by  atomic  oxygen 
slightly  favours  the  E-R  mechanism  since  the  adsorption  of  CO  with  the  presence  of  absorbed 
oxygen  is  very  weak  and  both  its  geometrical  and  electronic  structures  are  closer  to  that  of  a  free 
CO.  The  activation  energy  for  CO  oxidation  is  only  6.9  kcal/mol  by  molecular  oxygen  and  2.1 
kcal/mol  by  atomic  oxygen,  respectively.  So,  both  molecular  and  atomic  oxygen  are  very 
reactive  for  CO  oxidation  on  Ag(lOO).  The  formation  and  removal  of  CO3  on  silver  is  also 
discussed,  but  this  pathway  is  not  favoured  thermodynamically.  However,  the  low  adsorption 
energy  of  CO  on  Ag(lOO)  implies  a  low  sticking  probability  of  CO  on  silver  with  or  without  the 
presence  of  oxygen,  which  is  not  considered  in  the  current  study  as  well  as  the  effects  of 
temperature. 


E5:  References 

(1)  Liu,  L.;  Zhou,  F.;  Wang,  L.;  Qi,  X.;  Shi,  F.;  Deng,  Y.  Low -temperature  CO  oxidation  over 
supported  Pt,  Pd  catalysts:  Particular  role  of  FeOx  support  for  oxygen  supply  during 
reactions.  J.  Catal.  2010,  274,  1-10. 

(2)  Okumura,  M.  and  Haruta.  M.  Interplay  of  theoretical  calculations  and  experiments  for  a 
study  of  catalysis  by  gold.  Catal.  Today  2015,  259,  81-86. 

(3)  Daniel,  M.-C.;  Astruc,  D.  Gold  nanoparticles:  assembly,  supramolecular  chemistry, 
quantum-size-related  properties,  and  applications  toward  biology,  catalysis,  and 
nanotechnology.  Chem.  Rev.  2004, 104,  293-346. 

(4)  Min,  B.  K.;  Friend,  C.  M.  Heterogeneous  Gold-Based  Catalysis  for  Green  Chemistry:  Low- 
Temperature  CO  Oxidation  and  Propane  Oxidation.  Chem.  Rev.  2007, 107,  2709-2724. 

(5)  Luo,  J.;  Meng,  M.;  Li,  X.;  Zha,  Y.;  Hu,  T.;  Xie,  Y.;  Zhang,  J.  Mesoporous  Co304-Ce02  and 
Pd/Co304-Ce02  catalysts:  Synthesis,  characterization  and  mechanistic  study  of  their 
catalytic  properties  for  low -temperature  CO  oxidation.  J.  Catal.  2008,  254,  310-324. 

(6)  Zheng,  X.;  Zhang,  X.;  Wang,  X.;  Wang,  S.;  Wu,  S.  Low -Temperature  CO  Oxidation  over  a 
Ternary  Oxide  Catalyst  with  High  Resistance  to  Hydrocarbon  Inhibition.  Appl.  Catal.  A 
2005,  295,  142-149. 

(7)  Binder,  A.  J.;  Toops,  T.  J.;  Unocic,  R.  R.;  Parks  II,  J.  E.  and  Dai,  S.  Low -Temperature  CO 
Oxidation  over  a  Ternary  Oxide  Catalyst  with  High  Resistance  to  Hydrocarbon  Inhibition. 
Chem.  Int.  Ed.  2015,  54,  13263-3267. 

(8)  Ko,  E.-Y.;  Park,  E.  D.;  Lee,  H.  C.;  Lee,  D.  and  Kim,  S.  Supported  Pt-Co  Catalysts  for 
Selective  CO  Oxidation  in  a  Hydrogen-Rich  Stream.  Angew.  Chem.  Int.  Ed.  2007,  46,  734- 
737. 

(9)  Mu,  R.;  Fu,  Q.;  Xu,  H.;  Zhang,  H.;  Huang,  Y.;  Jiang,  Z.;  Zhang,  S.;  Tan,  D.  and  Bao,  X. 


52 


PROPOSAL  NUMBER:  62940-CH-REP,  FINAL  REPORT  2013-2017 


Synergetic  Effect  of  Surface  and  Subsurface  Ni  Species  at  Pt-Ni  Bimetallic  Catalysts  for  CO 
Oxidation.  J.  Am.  Chem.  Soc.  2011, 133,  1978-1986. 

(10)  Alayoglu,  S.;  Nilekar,  A.  U.;  Mavrikakis  ,  M.;  Eichhom,  B.  Ru-Pt  core-shell  nanoparticles 
for  preferential  oxidation  of  carbon  monoxide  in  hydrogen.  Nat.  Mater.  2008,  7,  333-338. 

(11)  Yamamoto,  T.  A.;  Nakagawa,  T.;  Seino,  S.;  Nitani,  H.  Bimetallic  nanoparticles  of  PtM  (M= 
Au,  Cu,  Ni)  supported  on  iron  oxide:  Radiolytic  synthesis  and  CO  oxidation  catalysis. 
Appl.Catal.  A  2010,  387,  195-202. 

(12)  Zhang,  H.;  Eiu,  X.;  Zhang,  N.;  Zheng,  J.;  Zheng,  Y.;  Ei,  Y.;  Zhong,  C.-J.;  Chen,  B.  H. 
Construction  of  ultrafine  and  stable  PtEe  nano-alloy  with  ultra-low  Pt  loading  for  complete 
removal  of  CO  in  PROX  at  room  temperature.  Appl.  Catal.  B  2016, 180,  237-245. 

(13)  Ma,  D.;  Tang,  Y.;  Yang,  G.;  Zeng,  J.;  He,  C.;  Eu,  Z.  CO  catalytic  oxidation  on  iron- 
embedded  monolayer  M0S2.  Appl.  Surf.  Sci.  2015,  328,  71-77. 

(14)  Ma,  D.  W.;  Ei,  T.;  Wang,  Q.;  Yang,  G.;  He,  C.;  Ma,  B.;  Eu,  Z.  Graphyne  as  a  promising 
substrate  for  the  noble-metal  single-atom  Catalysts.  Carbon  2015,  95,  756-765. 

(15)  Ei,  Y.;  Zhou,  Z.;  Yu,  G.;  Chen,  W.  and  Chen,  Z.  CO  Catalytic  Oxidation  on  Iron-Embedded 
Graphene:  Computational  Quest  for  Eow-Cost  Nanocatalysts.  J.  Phys.  Chem.  C  2010,  114, 
6250-6254. 

(16)  Song,  E.  H.;  Wen,  Z.;  Jiang,  Q.  CO  catalytic  oxidation  on  copper-embedded  graphene,  J. 
Phys.  Chem.  C2011, 115  (9),  3678-3683. 

(17)  Eiu,  X.;  Sui,  Y.;  Duan,  T.;  Meng,  C.;  Han,  Y.  CO  oxidation  catalyzed  by  Pt  embedded 
graphene:  a  first-principles  investigation,  Phys.  Chem.  Chem.  Phys.  2014,  16  (43),  23584- 
23593. 

(18)  Ei,  E.;  Zhao,  J.  and  Chen,  Z.  Ee-Anchored  Graphene  Oxide:  A  Eow-Cost  and  Easily 
Accessible  Catalyst  for  Eow-Temperature  CO  Oxidation.  J.  Phys.  Chem.  C  2012,  116, 
2507-2514. 

(19)  Mao,  K.;  Ei,  E.;  Zhang,  W.;  Pei,  Y.;  Zeng,  X.  C.;  Wu,  X.  et  al.  A  theoretical  study  of  single¬ 
atom  catalysis  of  CO  oxidation  using  Au  embedded  2D  h-BN  monolayer:  a  CO-promoted 
02  activation,  Sci.  Rep.  2014,  4,  1-7. 

(20)  Ein,  S.;  Ye,  X.;  Johnson,  R.  S.;  Guo,  H.  Eirst-principles  investigations  of  metal  (Cu,  Ag,  Au, 
Pt,  Rh,  Pd,  Ee,  Co,  and  Ir)  doped  hexagonal  boron  nitride  nanosheets:  stability  and  catalysis 
of  CO  oxidation,  J.  Phys.  Chem.  C2013, 117(33),  17319-17326. 

(21)  Eiu,  X.;  Duan,  T.;  Sui,  Y.;  Meng,  C.;  Han,  Y.  Copper  atoms  embedded  in  hexagonal  boron 
nitride  as  potential  catalysts  for  CO  oxidation:  a  first-principles  investigation.  RSC  Adv. 
2014,4(73),  38750-38760. 

(22)  Sinthika,  S.;  Kumar,  E.  M.;  Thapa,  R.  Doped  h-BN  monolayer  as  efficient  noble  metal-free 
catalysts  for  CO  oxidation:  the  role  of  dopant  and  water  in  activity  and  catalytic  de¬ 
poisoning,  J.  Mater  Chem.  A  2014,  2  (32),  12812-12820. 

(23)  Eiu,  X.;  Duan,  T.;  Meng,  C.;  Han,  Y.  Pt  atoms  stabilized  on  hexagonal  boron  nitride  as 
efficient  single-atom  catalysts  for  CO  oxidation:  a  first-principles  investigation.  RSC  Adv. 
2015,5(14),  10452-10459. 


53 


PROPOSAL  NUMBER:  62940-CH-REP,  FINAL  REPORT  2013-2017 


(24)  Narasimharao,  K.;  Al-Shehri,  A.;  Al-Thabaiti,  S.  Porous  Ag-Fe203nanocomposite  catalysts 
for  the  oxidation  of  carbonmonoxide.  Appl.  Catal.  A  2015,  505,  431-440. 

(25)  Afanasev,  D.  S.;  Yakovina,  O.  A.;  Kuznetsova,  N.  L;  Lisitsyn,  A.  S.  High  activity  in  CO 
oxidation  of  Ag  nanoparticles  supported  on  fumed  silica.  Catal.  Commun.  2012,  22,  43-47. 

(26)  Zhang,  X.  D.;  Qu,  Z.  P.;  Yu,  F.  L.;  Wang,  Y.  High-temperature  diffusion  induced  high 
activity  of  SBA-15  supported  Ag  particles  for  low  temperature  CO  oxidation  at  room 
temperature.  J.  Catal.  2013,  297,  264-271. 

(27)  Yousefi,  N.;  Pazouki,  M.;  Alizadeh,  M.  and  Hesari,  F.  A.  Biosynthesis  of  Ag/MWf-CNT 
Nanocomposites  Using  Aspergillus  fumigatus  as  the  CO  Oxidation  Catalyst.  Synthesis  and 
Reactivity  in  Inorganic,  Metal-Organic,  and  Nano-Metal  Chemistry  2016,  46,  464-470. 

(28)  Zhang,  X.  D.;  Qu,  Z.  P.;  Li,  X.  Y.;  Wen,  M.;  Quan,  X.;  Ma,  D.;  Wu,  J.  J.  Studies  of  silver 
species  for  low-temperature  CO  oxidation  on  Ag/Si02  catalysts.  Sep.  Purif.  Technol.  2010, 
72,  395-400. 

(29)  Tang,  D.;  Chen,  Z.;  Hu,  J.;  Sun,  G.;  Lu,  S.  and  Hu,  C.  CO  oxidation  catalyzed  by  silver 
nanoclusters:  mechanism  and  effects  of  charge.  Phys.  Chem.  Chem.  Phys.  2012,  14,  12829- 
12837. 

(30)  Chi,  Y.;  Zhao,  L.;  Lu,  X.;  An,  C.;  Guo,  W.;  Wu  C.  L.  Effect  of  alloying  on  the  stabilities 
and  catalytic  properties  of  Ag-Au  bimetallic  subnanoclusters:  a  theoretical  investigation.  J. 
Mater.  Sci.  2016,  51,  5046-5060. 

(31)  Negreiros,  F.  R.;  Sementa,  L.;  Barcaro,  G.;  Vajda,  S.;  Apra,  E.  and  Eortunelli,  A.  CO 
Oxidation  by  Subnanometer  Ag^Aua-^  Supported  Clusters  via  Density  Eunctional  Theory 
Simulations.  ACS  Catal.  2012,  2,  1860-1864. 

(32)  Wei,  Z.;  Ei,  D.;  Pang,  X.;  Ev,  C.  and  Wang,  G.  The  Mechanism  of  Eow-Temperature  CO 
Oxidation  on  IB  Group  Metals  and  Metal  Oxides.  ChemCatChem  2012,  4,  100-1 11. 

(33)  Gaussian  09,  Revision  B.Ol,  Erisch,  M.  J.  et  ah,  Gaussian,  Inc.,  Wallingford  CT,  2010. 

(34)  Perdew,  J.  P.;  Chevary,  J.  A.;  Vosko,  S.  H.;  Jackson,  K.  A.;  Pederson,  M.  R.;  Singh,  D.  J.; 
Eiolhais,  C.  Atoms,  molecules,  solids,  and  surfaces:  Applications  of  the  generalized  gradient 
approximation  for  exchange  and  correlation.  Phys  Rev  B  1992,  46,  6671. 

(35)  Perdew,  J.  P.;  Burke,  K.;  Wang,  Y.  Generalized  gradient  approximation  for  the  exchange- 
correlation  hole  of  a  many-electron  system.  Phys  Rev  B  1996,  54,  16533. 

(36)  Soo,  T.;  Qin,  C.  Adsorption  and  Dissociation  of  Carbon  Trioxide  on  Ag(lOO).  Int.  J. 
Quantum  Chem.  2010, 110,  946-952. 

(37)  Kuang,  X.;  Wang,  X.;  Eiu,  G.  A  density  functional  theory  study  on  the  Ag„H  (n  =  1-10) 
clusters.  Struct  Chem.  2011,  22,  517-524. 

(38)  Kuang,  X.;  Wang  X.;  Eiu,  G.  A  density  functional  study  on  the  Au^Ag  (n  =  1-12)  alloy 
clusters.  J.  Alloys  Compd.  2013,  570,  46-56. 

(39)  Kuang,  X.;  Wang  X.;  Eiu,  G.  The  adsorptions  of  silver-doped  small  gold  clusters  toward 
carbon  monoxide  molecule.  Struct  Chem.  2012,  23,  671-679. 

(40)  Stewart,  D.  M.;  Mavros,  M.  G.  and  Micha  D.  A.  Eight  Absorption  by  Crystalline  and 
Amorphous  Silicon  Quantum  Dots  with  Silver  Adsorbates  and  Dopants.  J.  Phys.  Chem.  C 


54 


PROPOSAL  NUMBER:  62940-CH-REP,  FINAL  REPORT  2013-2017 


2012,116,  23107-23112. 

(41)  Hay,  P.  J.;  Wadt,  W.  R.  Ab  initio  effective  core  potentials  for  molecular  calculations. 
Potentials  for  the  transition  metal  atoms  Sc  to  Hg.  J  Chem.  Phys.  1985,  82,  270. 

(42)  Wadt,  W.  R.;  Hay,  P.  J.  Ab  initio  effective  core  potentials  for  molecular  calculations. 
Potentials  for  main  group  elements  Na  to  Bi.  J  Chem.  Phys.  1985,  82,  284. 

(43)  Hay,  P.  J.;  Wadt,  W.  R.  Ab  initio  effective  core  potentials  for  molecular  calculations. 
Potentials  for  K  to  Au  including  the  outermost  core  orbitals.  J  Chem.  Phys.  1985,  82,  299. 

(44)  Rassalov,  V.  A.;  Pople,  J.  A.;  Ratner,  M.;  Windus,  T.  L.  6-31G*  basis  set  for  atoms  K 
through  Zn.  J  Chem.  Phys.  1998, 109,  1223. 

(45)  Hariharan,  P.  C.;  Pople,  J.  A.  The  influence  of  polarization  functions  on  molecular  orbital 
hydrogenation  energies.  Theor.  Chim.  Acta  1973,  28,  213. 

(46)  Francl,  M.  M.;  Pietro,W.  J.;  Hehre,W.  J.;  Binkley,  J.  S.;  Gordon,  M.  S.;  DeFrees,  D.  J.; 
Pople,  J.  A.  Self-consistent  molecular  orbital  methods.  XXIII.  A  polarization-type  basis  set 
for  second-row  elements.  J  Chem.  Phys.  1982,  77,  3654. 

(47)  Fukui,  K.  The  path  of  chemical -reactions  -  The  IRC  approach.  Acc.  Chem.  Res.,  1981,  14, 
363-68. 

(48)  Hratchian  H.  P.  and  Schlegel,  H.  B.  Using  Hessian  updating  to  increase  the  efficiency  of  a 
Hessian  based  predictor-corrector  reaction  path  following  method.  J.  Chem.  Theory  and 
Comput.,  2005, 1,  61-69. 

(49)  Vosko,  S.  H.;  Wilk,  L.;  Nusair,  M.  Accurate  spin-dependent  electron  liquid  correlation 
energies  for  local  spin  density  calculations:  a  critical  analysis.  Can.  J.  Phys.  1980,  58,  1200. 

(50)  Lee,  C.;  Yang,  W.;  Parr,  R.  G.  Development  of  the  Colle-Salvetti  correlation-energy 
formula  into  a  functional  of  the  electron  density.  Phys.  Rev.  B  1988,  37,  785. 

(51)  Becke,  A.  D.  Density-functional  exchange-energy  approximation  with  correct  asymptotic 
behavior.  Phys.  Rev.  A  1988,  38,  3098. 

(52)  Gajdos,  M.;  Eichler  A.  and  Hafner,  J.  CO  adsorption  on  close-packed  transition  and  noble 
metal  surfaces:  trends  from  ab  m/do  calculations.  J.  Phys.:  Condens.  Matter,  2004,  16, 
1141-1164. 

(53)  Hansen,  W.;  Bertolo,  M.;  Jacobi,  K.  Physisorption  of  CO  on  Ag(lll):  investigation  of 
the  monolayer  and  the  multilayer  through  HREELS,  ARUPS,  and  TDS.  Surf.  Sci.  1991,  253, 
1. 

(54)  McElhiney,  G.;  Papp,  H.;  Pritchard,  J.  The  adsorption  of  Xe  and  CO  on  Ag  (111).  Surf.  ScL, 
1976,  54,  617 

(55)  Backx,  C.;  de  Groot,  C.  P.  M.  and  Biloen,  P.  Adsorption  of  oxygen  on  Ag  (110)  studied  by 
high  resolution  EES  and  TPD.  Surf.  Sci.,  1981, 104,  300-317. 

(56)  Qin,  C.;  Sremaniak,  E.  S.  and  Whitten,  J.  E.  CO  Adsorption  on  Ag(lOO)  and  Ag/MgO(100). 
J.  Phys.  Chem.  B  2006, 110,  11272-11276. 

(57)  Engle,  T.;  Ertl,  G.  A  molecular  beam  investigation  of  the  catalytic  oxidation  of  CO  on 
Pd(lll).  7.  Chem.  Phys.  1978,  69,  1267. 


55 


PROPOSAL  NUMBER:  62940-CH-REP,  FINAL  REPORT  2013-2017 


(58)  Engle,  T.  A  molecular  beam  investigation  of  He,  CO,  and  O2  scattering  from  Pd(lll).  J. 
Chem.  Phys.  1978,  69,  373. 

(59)  Kuipers,  E.  W.;  Vardi,  A.;  Danon,  A.;  Amirav,  A.  Surface-molecule  proton  transfer:  A 
demonstration  of  the  Eley-Rideal  mechanism.  Phys.  Rev.  Lett.  1991,  66,  116-119. 

(60)  Rettner,  C.  T.  Dynamics  of  the  direct  reaction  of  hydrogen  atoms  adsorbed  on  Cu(lll) 
with  hydrogen  atoms  incident  from  the  gas  phase.  Phys.  Rev.  Lett.  1992,  69,  383-386. 

(61)  Stampfl,  C.;  Scheffler,  M.  Anomalous  Behavior  of  Ru  for  Catalytic  Oxidation:  A 
Theoretical  Study  of  the  Catalytic  Reaction  C0-i-l/202^C02.  Phys.  Rev.  Lett.  1997,  78, 
1500-1503. 


56 


PROPOSAL  NUMBER:  62940-CH-REP,  FINAL  REPORT  2013-2017 


Novel  MC/BZY  Proton  Conductor:  Materials  Development,  Device  Evaluation, 
and  Theoretical  Exploration  using  Cl  and  DFT  Methods 

Changyong  Qin,  Benedict  College,  Columbia,  SC  29204  (PI) 

Kevin  Huang,  University  of  South  Carolina,  Columbia,  SC  29208  (Co-PI) 

(Education  and  Community  Services  Part) 


PART  A:  STUDENT  MENTORING 

The  following  students  have  been  involved  in  the  research  activities  supported  by  the  current 
grant  in  the  year  of  2013-2017.  Six  of  the  seven  students  were  female  African  America  and  the 
work  funded  from  this  grant  had  served  as  their  requirements  for  senior  thesis/oral  defense.  One 
graduate  student  completed  a  Ph.D.  degree  with  partial  support  from  the  grant  in  the  Co-PI  group 
at  the  University  of  South  Carolina 

Kahla  Haines,  Chemistry  Major,  Benedict  College,  B.S.  Degree  2015 
Susan  Njoki,  Chemistry  Major,  Benedict  College,  B.S.  Degree  2016 
June  P.  Medina,  Chemistry  Major,  Benedict  College,  B.S.  Degree  2017 
Renee  Terrell,  Biology  Major,  Benedict  College,  B.S.  Degree  2017 
Shanice  Fezeu,  Biology  Major,  Benedict  College,  B.S.  Degree  2017 
Tamia  Brice,  Chemistry  Major,  Benedict  College,  B.S.  Degree  2019  (expected) 

Maria  Melville,  Biology  Major,  Benedict  College,  B.S.  Degree  2019  (expected) 

Xiaolei  Xiong,  Mechanical  Engineering,  PhD  in  2017,  University  of  South  Carolina 
Dr.  Xueling  Lei,  Postdoc  Research  Associate 
Dr.  Jingjing  Tong,  Postdoc  Research  Associate 
PART  B:  INFRASTRUCTURE  IMPROVEMENT 

The  Undergraduate  Teaching  and  Research  Laboratory  for  Advanced  Materials  (UTRAM)  has 
been  developed  with  increased  research  capacity  and  is  now  the  Center  for  Advanced  Materials 
(CAM)  at  Benedict  College.  It  serves  a  core  facility  for  faculty/student  research  and  also  for  the 
STEM  education  and  community  services.  In  particular,  follow  instruments  have  been  purchased 
and  added  to  CAM  through  the  current  grant  supports. 

Vacuum  drying  oven  (QTN:  1) 

VERSASTA3-200  Advanced  DC  Voltammetry  System  (QTN:  1) 

Solartron  1470  Luel  Cell  Test  System  (QTN:  1) 
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High  temperature  sample  stage  for  Raman  Spectroscopy  (QTN:  1) 

64-Core  Intel  Computing  Server  (QTN:  1) 

AMD  Opteron  Quad-Socket  Computing  System  (QTN:  4) 

Rotating  Disc  Electrode  System  (QTN:  1) 

Planetary  Ball  Mill  (QTN:  1) 

Rotary  Ball  Mill  (QTN:  1) 

PART  C:  PUBLICATIONS  AND  PRESENTATIONS 

The  project  funded  by  the  current  grant  has  led  to  several  conference  and  meeting  presentations, 
and  also  six  papers  published  in  peer-reviewed  journals  and  two  manuscripts  being  submitted. 
Below  is  the  complete  list. 

CO  Oxidation  by  Molten  Carbonate:  A  DFT  Study 
S.  Njoki  and  C.  Qin,  WoPhyS’15,  Lincoln,  Nebraska 

Proton  Transfer  in  Molten  Carbonate:  Mechanism  and  Kinetics  from  DFT  Study 
X.  Lei  and  C.  Qin,  APS  annual  spring  meeting,  Washington,  DC,  March,  2016 

Molten  Carbonate  as  Non-metal  Catalyst  for  CO  Oxidation 
J.  Tong  and  C.  Qin  manuscript  in  submission. 

Kinetics  and  Mechanism  of  Proton  Transfer  in  Molten  Lithium  Carbonate:  Insights  from 
Static  and  Dynamic  DFT  Studies 

X.  Lei,  K.  Huang,  C.  Qin,  Scientific  Reports,  DOT  I0.I038/s4I598-0I7-07726-3. 

DFT  Study  of  CO  Oxidation  by  Molten  Carbonate 
X.  Lei,  K.  Huang,  C.  Qin,  manuscript  in  preparation. 

DFT  Study  of  Oxygen  Migration  in  Molten  Carbonate 

X.  Lei,  K.  Haines,  K.  Huang,  C.  Qin,  J.  of  Power  Sources,  2016,  305,  161-166. 

Enhanced  Interfacial  Proton  Migration  on  BaZr(Y)03  by  Molten  Carbonate:  A  First 
Principles  Study 

X.  Lei,  K.  Huang,  C.  Qin,  Solid  State  Ionics  ,2016,  289,  48-54. 

DFT  Study  of  Oxygen  Dissociation  in  Molten  Carbonate 

X.  Lei,  K.  Haines,  K.  Huang,  C.  Qin,  J.  of  Phys.  Chem.  A,  2015,  119,  8806-8812. 

Energetics  of  Proton  Transfer  in  Alkali  Carbonates:  A  First  Principles  Calculation 
X.  Lei,  C.  Qin,  K.  Huang,  RSC  Advances,  2015,  5,  56205-56209. 

Synergetic  Proton  Conduction  in  BZY-Carbonate  Composite  Electrolyte  for  IT-SOFC 
X.  Xiong,  X.  Lei,  C.  Zhang,  J.  Wang,  K.  Huang,  Solid  State  Ionics  ,  2015,  279,  66-71. 

PART  C:  OUTREACH  ACTIVITIES 

The  grant  has  supported  the  Xtreme  Technology  Event  for  high  school  students  at  Benedict 
College.  Most  of  them  are  from  African  American  community.  Three  African  American 
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middle/high  school  students  have  also  participated  a  one-week  fuel  cell  program  at  the  Co-PI 
laboratory  at  The  University  of  South  Carolina. 

The  grant  has  supported  the  Xtreme  Technology  Event  for  high  school  students  at  Benedict 
College.  Most  of  them  are  from  African  American  community.  The  grant  also  allowed  the  STEM 
students  to  participate  the  annual  Harambee  festival  on  the  campus  of  Benedict  College,  which 
services  the  local  African  American  community  and  increases  their  awareness  of  high 
technologies. 


